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ABSTRACT 

The  goal  of  this  research  was  to  develop  analysis  and  optimal  design  procedures  for 
planar  as  well  as  spatial  mehcanisms  that  arc  frequently  used  in  space  structures.  A  nonlinear 
finite  element  procedure  which  was  developed  originally  for  planar  mechanisms  during  the  initial 
stages  of  this  research,  has  been  modified  considerably  to  handle  complex  mechanisms  with  sliding 
masses  and  mechanisms  operating  at  relatively  high  speeds.  The  analysis  takes  into  account  the 
effects  of  geometric  and  material  nonlinearities,  vibrational  effects  and  coupling  of 
deformations.  Numerical  results  have  been  reported  for  certain  mechanism  examples.  The  effects 
of  nonlinearities  have  been  found  to  be  significant  on  the  dynamic  behavior  of  mechanisms. 
Considerable  progress  has  been  made  in  developing  a  nonlinear  finite  element  procedure  for 
three-dimensional  mechanisms.  Numerical  results  obtained  for  some  example  problems  indicate  the 
validity  of  the  current  three-dimensional  formulation.  A  new  optimization  algorithm  has  also 
been  developed  based  on  the  Gauss  method  to  handle  various  types  of  nonlinear  constraints  with 
the  goal  of  reducing  the  number  of  analyses  required  to  obtain  an  optimal  design.  Complete 
details  of  the  nonlinear  finite  element  procedures  as  well  as  the  optimization  technique  are 
available  in  published  papers,  copies  of  which  are  included  here  in  the  Appendix.  Because  of  the 
complex  nature  of  the  nonlinear  analysis,  which  had  to  be  repeated  many  times  during  the 
optimization  process,  considerable  amount  of  computer  time  was  needed  for  this  research.  To  help 
overcome  these  computational  difficulties,  DoD  and  NSF  provided  funds  through  two  separate  grants 
to  purchase  a  research  computer  and  other  associated  equipment  as  well  as  access  to  a 
supercomputer. 
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SUMMARY 


A  nonlinear  finite  element  analysis  procedure  which  was  developed  for  planar  mechanisms 
during  the  first  phase  of  this  research  has  undergone  considerable  modification  to  handle  complex 
mechanisms  with  sliding  masses  and  mechanisms  operating  at  relatively  high  speeds.  A  suitable 
nonlinear  finite  element  analysis  procedure  has  also  been  developed  for  three-dimensional 
mechanisms.  In  both  cases,  the  analysis  takes  into  account  the  effects  of  geometric  and  material 
nonlinearities,  vibrational  effects,  and  coupling  of  deformations.  Numerical  results  have  been 
reported  for  certain  mechanism  examples.  These  results  indicate  significant  influences  of  the 
geometric  as  well  as  the  material  nonlinearity  effects  on  the  dynamic  behavior  of  mechanisms.  In  the 
optimal  design  area,  a  new  algorithm  has  been  developed  for  finding  the  minimum  of  a  sum-of-squares 
objective  function  subject  to  general  nonlinear  constraints.  The  solution  of  some  selected  examples 
indicate  good  results  in  terms  of  the  total  number  of  objective  function  evaluations  to  obtain  an 
optimal  design.  Complete  details  of  these  investigations  as  well  as  those  of  the  nonlinear  finite 
element  analysis  are  included  in  the  Appendix.  To  meet  the  extraordinary  computational  needs  of  this 
project,  a  separate  VAX  1 1/785  computer  and  peripheral  equipment  were  made  available  through  a  DoD 
research  grant.  The  National  Science  Foundation  also  provided  funds  for  some  additional  equipment  as 
well  as  computational  time  on  a  supercomputer.  The  current  research  will  continue  with  the  help  of  a 
recently  obtained  National  Science  Foundation  grant  (Grant  No.  INT-8616036)  which  involves  the 
development  of  a  joint  clearance  model  for  mechanical  mechanisms  and  its  inclusion  in  the  vibrational 
analysis  of  linkages  undergoing  large  deformations  due  to  high  speeds  and  large  loads. 

RESEARCH  OBJECTIVES 

The  objective  of  this  research  was  to  develop  a  nonlinear  finite  element  dynamic  analysis 
procedure  for  planar  as  well  as  spatial  mechanisms  that  are  frequently  used  in  space  structures. 

Included  in  the  nonlinear  analysis  are  the  effects  of  curvature-displacement  nonlinearity, 
nonlinearity  due  to  extension  or  streiching  (both  caused  by  large  defonmations),  material 
nonlinearity  as  well  as  combinations  of  these.  In  addition  to  the  nonlinear  analysis,  an  efficient 
optimal  design  method  was  to  be  developed  to  handle  objective  functions  composed  of  combinatioas  of 


rigid  body  and  deformation  displacements  involving  geometric  design  variables  as  well  as 
cross-sectional  sizes  of  the  members  of  the  mechanism  subject  to  limitations  on  stresses  and 
deformations.  Thus  the  proposed  research  involves  the  following. 

a)  develop  a  nonlinear  finite  element  procedure  for  dynamic  planar  mechanisms. 

b)  develop  an  efficient  optimization  method  involving  a  small  number  of  analyses  for  mechanism  design 

problems. 

c)  extend  the  nonlinear  analysis  procedure  developed  for  two-dimensional  mechanisms  to 
three-dimensional  mechanisms. 

d)  apply  the  optimization  technique  developed  to  simple  mechanism  problem. 

e)  combine  the  nonlinear  analysis  procedure  with  the  optimization  technique  to  design  complex 

three-dimensional  mechanisms,  robots  and  mechanical  manipulators. 

SIGNIFICANT  ACCOMPLISHMENTS 

A  nonlinear  finite  element  analysis  procedure  has  been  developed  for  planar  mechanisms  to 
handle  geometric  and  material  nonlinearities  (see  the  publications  list  and  Appendix  for  details). 

The  results  of  several  example  mechanisms  clearly  indicate  the  need  to  include  these  types  of 
nonlinearities  in  the  dynamic  analysis.  The  difficulties  in  extending  this  approach  to  complex 
planar  mechanisms  with  sliding  masses  and  mechanisms  operating  at  relatively  high  speeds  have  been 
overcome  by  using  a  modified  finite -element  formulation  to  handle  such  complex  cases.  Some  example 
problems  using  this  new  formulation  have  been  considered  in  a  paper  entitled,  "Vibrational  Analysis 
of  Mechanisms  with  Geometric  and  Material  Nonlinearities"  (with  E.  Kear,  M.  Sathyamoorthy  and  K.D. 
Willmert  as  authors)  presented  at  the  Joint  AFOSR-SES  Meeting  held  at  the  State  University  of  New 
York  at  Buffalo  in  August  1986.  Similar  results  have  also  been  presented  in  a  very  recent  paper 
entitled,  "Effect  of  Geometric  and  Material  Nonlincaritics  on  Vibration  of  Planar  Mechanisms"  by  E.B. 

Kear,  M.  Sathyamoorthy  and  K.D.  Willmert,  presented  at  the  ARO/AFOSR  Meeting  on  Nonlinear  Vibrations, 
Stability  and  Dynamics  of  Structures  and  Mechanisms  held  at  VPI  &  SU,  Blacksburg,  Virginia,  March 
1987.  Considerable  progress  has  also  been  made  in  developing  a  nonlinear  finite  clement  procedure 
for  three-dimensional  mechanisms.  Some  example  problems  have  been  solved  and  the  results  indicate 


the  validity  of  the  three-dimensional  formulation.  The  general  three-dimensional  formulation  is 
given  in  a  paper  entitled,  "Finite-Element  Nonlinear  Analysis  of  Three-Dimensional  Mechanisms."  by 
M.  El-Sawy,  K.D,  Willmert  and  M.  Sathyamoorthy  included  in  the  Appendix  of  this  report. 

In  the  optimization  area,  a  very  efficient  optimality  criterion  technique  called  the  Gauss 
Constrained  Method  has  been  developed  to  solve  optimal  design  problems  with  objective  functions  which 
are  the  sum  of  squared  quantities  with  general  nonlinear  constraints.  The  technique  has  the 
advantage  of  reducing  the  number  of  analyses  required  to  obtain  an  optimal  design,  thereby 
significantly  reducing  the  computational  time.  This  method  is  described  in  a  paper  entitled,  "The 
Development  and  Application  of  Gauss’  Nonlinearly  Constrained  Optimization  Method"  (with  D.R.  Boston, 
K.D.  Willmert  and  M.  Sathyamoorthy  as  authors)  published  in  the  Journal  of  Computer  Methods  in 
Applied  Mechanics  and  Engineering.  In  a  paper  entitled,  "The  Gauss  Optimization  Method  for  Problems 
with  General  Nonlinear  Constraints"  by  T.E.  Potter,  K.D.  Willmert  and  M.  Sathyamoorthy  (presented  at 
the  22nd  Annual  Meeting  of  the  Society  of  Engineering  Science  held  at  the  Pennsylvania  State 
University,  University  Park,  Pennsylvania,  October  1985),  a  new  algorithm  has  been  developed  for 
finding  the  minimum  of  a  sum-of-squares  objective  function  subject  to  general  nonlinear  constraints. 

The  solution  of  examples  indicate  good  results  in  terms  of  the  total  number  of  objective  function 
evaluations  required  by  the  algorithm  to  obtain  an  optimal  design.  The  optimization  techniques 
developed  in  this  research  as  extensions  of  the  Gauss  method  to  handle  various  types  of  constraints, 
reduce  the  number  of  analyses  required  to  obtain  an  optimal  design.  The  method  is  now  being  used  to 
solve  additional  example  problems  including  various  mechanisms. 

Because  of  the  very  complex  nonlinear  analysis  required,  which  must  be  repeated  many  times 
during  the  optimization  process,  a  considerable  amount  of  computer  time  was  needed  for  this 
research.  To  meet  these  needs,  a  proposal  entitled  "Laboratory  for  Graphical  Analysis  of  Nonlinear 
Deformations  in  Dynamic  Structural-Mechanical  Systems"  was  submitted  to  DoD  under  the  DoD  - 
University  Instrumentation  Program  to  purchase  a  separate  research  computer  for  this  project.  This 
resulted  in  a  grant  (No.  AFOSR-85-0103)  of  $101 ,567.  Although  the  original  proposal  called  for  the 
purchase  of  a  VAX  1 1/730,  a  very  careful  and  thorough  search  for  the  best  computer  (with  the 
available  funds)  resulted  in  the  purchase  of  a  much  larger  and  faster  VAX  1 1/785.  Digital  Equipment 


Corporation  offered  a  sizable  reduction  in  cost  of  its  VAX  1 1/785  computer  under  the  DEC  Educational 
Discount  Program.  Because  of  this  reduction  and  because  of  additional  cost  sharing  by  Clarkson 
University’s  School  of  Engineering,  it  was  possible  to  purchase  the  VAX  1 1/785  at  no  additional  cost 
to  DoD.  The  computer  equipment  purchased  through  DoD-URIP  Grant  and  other  sources  includes 


I.  VAX  1 1/785  COMPUTER  HARDWARE 

1.  VAX  1 1/780  Packaged  System  Including: 

(A)  VAX  11/780  CPU 

(B)  2-Mbytes  ECC  MOS  (64-K  chip)  Memory 

(C)  H9652  UNIBUS  Expansion  Cabinet  with  BA  11-K  and  DD11-DK 

(D)  VAX/VMS  License  and  Warranty 

2.  TU80  9  Track  Streaming  Tape  Drive  with  Cabinet 

3.  RUA81  456  Mbyte  Fixed  Disk 

4.  2-Mbytes  of  Additional  Memory 

5.  FP780  Floating  Point  Accelerator 

6.  Two  DMF32-LP  Communication  Interfaces 

7.  780  to  785  Upgrade  Kit 

8.  25  ft.  RS  232  Sync  Cable 

9.  Two  300/1200/2400  Baud  Telephone  Modems 

10.  Installation 

1 1 .  Insurance  and  Transportation 

12.  Miscellaneous  -  Installation  of  Power,  Phone 

Lines,  Terminal  Lines,  etc. 

Total  Computer  Hardware  Cost 

II.  COMPUTER  DISPLAY  TERMINAL 

1.  Tektronix  M41 15B  Computer  Display  Terminal 

2.  Option  N 1 :  Warranty-Plus 

3.  Option  2B:  Additional  512  Kbytes  RAM 

4.  Option  23:  Additional  Four  Planes  of  Display  Memory 


$102,750 


$8,800 

$19,600 

$8,100 

$8,960 

$5,250 

$80,000 


$1,060 


$1,817 


$1,756 


$238,188 


$19,950 


$1,025 


$4,600 


$6,000 


«  a.  x  _  *  -N  f  «.%  .  ■  «  "*  a.  »  _%  _  »  a*  *  a."*  -  *  k  »  k"*  a  II  a  V  .V  JV  k.11  a.%  a.*  -*S»  -  ■  -  "k  *.">  • 


$500 


5.  Option  09:  4695  Color  Copier  Interface 

6.  Option  42:  Single  Flexible  Disk  $1,700 

7.  Display  Stand  $750 

8.  Software  Package  $1,000 

9.  4695  Color  Graphic  Copier  $1,595 

10.  Option  42:  Warranty-Plus  $430 

11.  4926  10  Mbyte  Hard  Disk  $4,200 

12.  Option  Nl:  Warranty-Plus  $210 

13.  Shipping  $371 

Total  Display  Terminal  Cost  $42,33 1 

III.  SOFTWARE  FOR  VAX  1 1/785  COMPUTER 

1 .  VMS  Operating  System  N/C 

2.  FORTRAN  License  $5,170 

3.  DECNET  Communication  Software  $2,950 

4.  IGL  Graphics  Software  $2,677 

5.  PSI  Access  Software  $1,850 

Total  Software  Cost  $12,647 

IV.  TOTAL  HARDWARE  &  SOFTWARE  COSTS  $293,166 


The  total  value  of  the  equipment  and  software  is  $293,166.  Discounts  and  contributions  from 
Digital  Equipment  Corporation,  Tektronix,  Clarkson  University’s  College  of  Engineering  and  the 
Department  of  Mechanical  and  Industrial  Engineering  total  $191,599.  Thus,  the  total  cost  of  the 
hardware  and  software  to  DoD  remained  at  $101 ,567  as  originally  proposed.  It  should  be  noted  that 
the  capabilities  of  the  VAX  1 1/785  system,  including  the  Tektronix  4115  graphic  terminal,  arc 
enormous  compared  to  the  originally  proposed  VAX  1 1/730.  The  VAX  1 1/785  system  is  five  times  faster 
than  the  VAX  1 1/730,  has  4  Mbytes  of  memory  (compared  to  only  1  Mbyte  of  memory  for  VAX  1 1/730),  456 
Mbyte  of  disk  space  (compared  to  121  Mbyte  of  disk  space)  a  total  of  16  terminal  lines,  and  a  9  track 
streaming  tape  drive  (no  tape  drive  was  included  in  the  original  VAX  1 1/730  system). 

Hardware  and  software  were  also  purchased  to  tic  the  VAX  1  1/785  into  Clarkson’s  campus-wide 
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computer  network.  The  physical  link  is  through  the  University’s  VAX  1 1/780,  but  this  is  tied  to  the 
other  computers  on  campus,  which  is  linked  to  other  universities  though  BITNET.  This  tie  in  of  the 
VAX  1 1/785  allows  the  users  of  this  research  computer  access  to  many  of  the  other  facilities  of  the 
university,  such  as  high  speed  printers,  digital  plotters,  laser  printers,  etc.  It  also  allows 
researchers  with  terminal  connected  to  the  other  computers  on  campus  to  sign  on  to  the  VAX  1 1/785  as 
though  they  were  directly  connected. 

The  Tektronix  41 15B  computer  display  terminal,  which  is  connected  to  the  VAX  1 1/785  computer, 
has  recently  been  expanded  to  improve  its  capabilities.  Both  a  3-dimensional  wire  frame  and  a  shaded 
surface  option  have  been  added.  These  options  allow  the  terminal  to  locally  manipulate  3-dimensional 
objects,  such  as  rotating  them  in  3-dimensional  space,  removing  hidden  lines,  drawing  shaded 
surfaces,  etc.  These  expansions  result  in  this  terminal  being  equivalent  to  a  Tektronix  4129 
terminal,  which  is  the  most  recent  high-end  terminal  introduced  by  Tektronix.  The  total  cost  of 
these  options  was  $16, 475,  which  was  made  possible  through  contributions  from  Gleason  Foundation, 
Proctor  and  Gamble,  Tektronix,  Coming  Glass  Works  as  well  as  the  University’s  School  of 
Engineering. 

A  grant  from  the  National  Science  Foundation  (Grant  No.  DMC-8500627),  with  M.  Sathyamoorthy 
and  K.D.  Willmert  as  principal  investigator,  included  funds  totaling  $10,715  for  the  further 
expansion  of  the  graphic  facilities.  A  Tektronix  4692  color  graphics  copier,  a  Tektronix  4107  low 
resolution  graphic  terminal  and  two  Z-200  personal  computers  were  purchased  through  this  NSF  grant. 
These  purchases  complement  the  high  resolution  Tektronix  4115B  terminal  obtained  through  the 
DoD-University  Research  lastmmentation  Program.  In  addition  to  these  equipment  funds,  this  NSF 
grant  provided,  as  part  of  its  Cooperative  Program  on  the  Use  of  Supercomputers,  twenty-five  hours  of 
CPU  time  on  a  Cray  X-MP  supercomputer.  The  program  development  and  trial  runs  were  done  on  the  DoD 
funded  in  house  VAX  1  1/785  research  computer  with  final  runs  made  at  the  supercomputer  located  at  the 
University  of  Illinois. 

A  summary  of  funding  sources  including  the  1984-85  DoD-URIP  Grant  to  purchase  the  VAX  1 1/785 
research  computer  and  other  associated  equipment  (including  upgrades)  is  given  below: 

DoD  -  University  Research  Instrumentation  Program  $101,567 
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Digital  Equipment  Corporation  Contribution 
Tektronix  Discount  and  Contribution 


$120,852 


B 

u 


National  Science  Foundation 

Clarkson  University’s  College  of  Engineering  Contribution 

Department  of  Mechanical  &  Industrial  Engineering 

Clarkson’s  Educational  Resource  Center 

Gleason  Foundation 

Proctor  and  Gamble 

Coming  Glass  Works 

TOTAL 


$11,110 

$10,715 

$56,999 

$3,176 

$3,300 

$7,000 

$5,300 

$337 

$320,356 


As  a  result  of  contributions  from  all  of  these  sources,  the  total  value  of  the  equipment  within  this 
laboratory  exceeds  $300,000  for  an  investment  of  only  slightly  over  $100,000  from  DoD. 

A  recent  grant  from  the  National  Science  Foundation  (Grant  No.  1NT-8616036),  with  K.D. 
Willmert  and  M.  Sathvamoorthy  as  principal  investigator,  will  help  accomplish  all  the  remaining  goals 
of  the  current  AFOSR  research.  The  funded  NSF  research  involves  a  cooperative  effort  between  the 
principal  investigators  and  a  collaborator  at  the  Korea  Advanced  Institute  of  Science  and  Technology 
in  the  Republic  of  Korea.  The  particular  project  will  require  the  development  of  a  joint  clearance 
model  for  mechanical  mechanisms  and  its  inclusion  in  the  vibrational  analysis  of  linkages  undergoing 
large  deformations  due  to  high  speeds  and  large  loads.  Also  included  is  the  study  of  the  optimal 
design  of  counter  weights  to  reduce,  as  much  as  possible,  the  joint  forces.  The  duration  of  the  NSF 
research  will  be  until  December  1989. 
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4.  INTEGRATED  DESIGN  OF  STRUCTURE  AND  CONTROL 


by  Vipperla  B.  Venkayya 


I  ntroduction 

The  panel  identified  several  presentations  relevant  to  the  subject 

of  integrated  design.  They  included  the  following: 

1.  Ken  Willmert  of  Clarkson  University  on  the  optimization  of  flexible 

mechanisms ; 

2.  Mohan  Aswani  and  G.  T.  Tseng  of  The  Aerospace  Corporation  on 

continuum  modeling  of  the  plant  as  one  means  of  simplifying 
structure-control  optimization  problems; 

3.  Manohar  Kamat  of  VPI  on  the  issues  of  plant  nonlinearities 

(geometric  and  material),  proposing  that  they  be  considered  in 
control  system  design; 

4.  Moktar  Salama  of  JPL  on  structure-control  optimization  with  a  single 
performance  index  consisting  of  the  structural  mass  and  the  total 
control  input; 

5.  Dale  Berry  of  Purdue  University  on  continuum  modeling  as  a  means 
of  reducing  plant  dimensionality; 

6.  K.  C.  Park  of  Lockheed  Palo  Alto  Research  Lab  on  the  use  of 

transient  energy  density  profiles  to  achieve  optimum  disturbance 

dissipation  and  control; 

7.  John  Junkins  of  VPI  on  some  recent  results  of  eigenvalue  placement 
via  structural  parameter  optimization. 

The  material  of  these  presentations  was,  to  a  certain  extent,  the  basis 

for  the  panel  discussions. 

The  panel  proposed  the  following  five  topics  for  discussion: 


49 


2 


Volume  1 
ENGINEERING 
MECHANICS 
IN  CIVIL  ENGINEERING 


Proceedings  ot  the 

Fittn  Engineering  Mechanics  Division 

Specialty  Conference 

Sponsored  by  the 

Engineering  Mechanics  Division 

of  the  American  Society  ot  Civil  Engineers 

University  ot  Wyoming 
Laramie  Wyoming  82071 
August  1-3.  1984 

in  cooperation  with 

Tne  Department  of  Civil  Engineering  and 
The  Division  of  Conferences  and  Institutes 
University  of  Wyoming 

Edited  by 

A.  P  Boresi  and  K.  P  Chong 


Tv  the 

A'TtP'i r.iP  N.oety  :jt  I  :■»"  £'3g  '>er’ 

•45  f-  jst  4  St'eef 

VjU* 


,v.v 


s  1  t 


optima:,  design 

5-  DYNAMIC  MECHANICAL  SYSTEMS 
•• SDF.RG01SG  URGE  DEFORMATIONS 

K.  D.  Wil Inert 
M.  Sathyamoor thy 

Mechanical  and  Industrial  Engineering  Department 
Clarkson  University 
Potsdam,  NY  13676 


Introduct ion 

Tr.e  general  purpose  of  most  mechanical  devices 
is  to  produce  accurate  two  or  three-dimensional 
movement  of  objects.  However,  tor  large  external 
loads  and/or  nigh  speed  operation,  sufficient  forces, 
both  inertial  and  externally  applied,  are  produced  to 
cause  their  members,  joint*,  and  support  structure  to 
detorm.  Ihi6  results  in  a  loss  of  accuracy  of  the 
device.  To  solve  this  problem,  the  current  design 
procedure  is  to  increase  the  stiffness  of  the  members 
by  modifying  their  cross-sect lona 1  sizes  (commonly  by 
increasing  the  areas)  or  changing  the  material  used 
to  reduce  the  deformations  and  stresses.  The  result 
is  a  more  massive  device  which  is  difficult  to 
. jntrol,  requiring  more  power  to  drive  and 
interacting  with  its  structural  support  to  an  even 
greater  extent . 

Optimization  methods  have  been  applied  to 
determine  the  member  sizes  by  minimizing  the  weight 
subject  to  limitations  on  deformations  and  stresses. 
This  can  reduce  the  overall  size  of  Che  device,  but 
result  is  still  not  the  best  design.  The 
: : : :  i:u I t y  lies  m  the  basic  separation  of  variables 
' ^  design  process.  The  device  configuration 
the  number  and  type  of  members  and  joints)  and 
::e  member  lengths  are  designed  fir6t  to  produce  the 
desirej  motion,  assuming  the  members  are  rigid. 
Trer,  independently.  the  cross-sect  ions l  sizes  are 
determined  so  that  the  actual  motion  deviates  as 
litt.p  as  possible  from  the  rigid  body  motion.  In 
t "  i  s  pr  edure,  the  deformations  of  the  members  are 
:or.sidered  to  na-.e  a  detrimental  effect  on  the 
jvprali  motion.  The  goal,  therefore,  is  to  restrict 
-r  reduce  them.  owever ,  these  deformations  can 
i:tually  be  used  to  mprove  the  motion  of  the  device 
it  tr.e  geometric  design  variables  (member  lengths, 
•'  t  and  c  ros  s- sec  c  iona  l  sizes  of  the  members  are 
mb  i  tied  and  treated  as  a  single  set  of  design 
variables.  Thus,  the  combined  design  problem  is  to 
st-le:t  ail  of  the  design  variables  so  that  the  actual 
roti  n  > f  the  device,  which  is  a  combination  of  the 
rigid  body  motion  plus  the  deformations,  is  as  close 
as  pa  ssible  to  the  desired  motion.  The  usual  rigid 
“j ■  -  iy  constraints  exist  on  thi6  problem,  such  as 
limitations  on  the  moveability  of  the  device, 

.  1 1  ;  n  s  of  the  support  points,  lengths  of  the 

t.  r  s  ,  etc.  But  also  present  are  deformation  and 

stress  .or.it  rainti ,  and  natural  frequency  limitations 
ass  . v  isted  with  the  flexibility  of  the  members. 

er.di  ng  :r,  t  tv  type  of  motion  desired,  speed 

■  :  :  *•  ration  of  the  device,  and  restrictions  on  the 

>,‘i  ur ,  the  ni'D.jt  i  m  A  this  design  problem  will  bp  a 


mechanism  in  which  the  deformations  of  the  members 
serve  to  aid  the  mechanism  in  producing  the  desired 
motion.  Thus  a  much  wider  range  of  motions  is 
possible  than  could  be  obtained  from  a  rigid  link 
device.  The  optimal  design,  however,  is  likely  to 
contain  several  members  which  are  very  flexible  (any 
members  that  are  required  to  be  rigid  can  be  designed 
using  appropriate  constraints).  Because  of  these 
large  deformations,  a  general  nonlinear  vibrational 
deformation  and  stress  analysis  is  required.  The 
types  of  non  1  inesr i t les  include: 

1.  Nonlinear  curvature-displacement  relations 

2.  Extension  or  stretching  of  the  neutral  surface 

3.  Material  nonlinearities  (nonlinear  stress- 
strain) 

4.  Effects  of  transverse  shear  and  rotatory 
inertia  due  to  realistically  proportioned  members 

5.  Joint  clearances 

6.  Variable  nature  of  the  cross-section  .if  the 
members 

7.  Interaction  of  the  mechanical  device  with  the 
support  structure 

8.  Coupling  of  the  degrees-of-f reedom  due  to  large 
de format  ions . 

A  typical  analysis  which  takes  into  account  all  of 
these  nonl  inesr it le s  is  iterative  in  nature  requiring 
several  iterations  to  obtain  a  solution. 

Because  of  the  complex  nonlinear  analysis 
required.  which  must  be  done  many  times  during  the 
optimization  phaae,  considerable  computational  time 
is  needed  for  this  design  procedure.  Experience  has 
shown  that  four  to  eight  hours  of  computer  time  is 
common  to  perform  just  one  analysis  of  a  single 
planar  four-bar  or  slider  crank  mechanism  on  an  IBM 
4341  computer.  For  more  complex  planar  mechanical 
devices  consisting  of  additional  links,  gears,  cams 
and  other  elements,  considerably  longer  times  would 
be  expected.  For  three-dimensional  mechanisms,  a6 
exist  in  many  applications  of  mechanical 
man lpu l a t or s ,  automotive  suspensions,  etc.,  enormous 
computational  times  would  be  invo 1 ved--aga in  to 
perform  just  one  analysis. 

The  mathematical  optimization  techniques 
currently  available  for  solving  optimal  design 
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problem#  all  require  several  iteration#  to  obtain  the 
test  den<n.  Some  method#  involve  a  large  number  of 
iteration#,  with  each  iteration  requiring  numerous 
analyse#  'obtaining  numerical  values  for  the 
•‘k  ;ect  ive  function  and  constraints  for  particular 
values  of  the  design  variables).  The#e  methods  are 
adequate  if  the  deiign  problem  is  small.  since 
. ooput at  ions l  times  are  relatively  insignificant. 
However  fcr  large  design  problem#,  or  ones  in  which  a 
complex  analysis  vs  required.  it  is  extremely 
important  that  the  optimization  technique  be 
efficient.  particularly  in  term#  of  the  number  of 
analyses  required.  Even  then,  the  solution  of  these 
problems  i#  impossible  on  most  computers  currently  ar 
universities.  The  availability  of  a  super  computer 
to  accomplish  the  optimal  design  of  large  mechanical 
svatems  is  absolutely  necessary. 

Reaearch  Problems 

The  particular  research  problems  in  the  area  of 
optimal  design  of  dynamic  mechanical  systems  which 
require  the  use  of  super  computers  are: 

•  Nonlinear  Ana l vs  is  of  Meehan  1  sms 

Objective:  Develop  s  large  amplitude  nonlinear 

deformation,  stress  and  frequency  analysis  of  two  and 
three-dimensional  mechanical  systems  undergoing  high 
speed  motion  subjected  to  large  external  loads. 

Include  : 

*  Geometric  and  material  nonl inear  it  le s 

*  The  effects  of  transverse  shear  and  rotatory 
inert  is 

*  Joint  clearances 

*  Interaction  of  the  mechanical  devices  with 
the  support  structure 

*  Coupling  of  degreea-of-f reedom  due  to  large 
de f ormat  ions 

Supercomputer  needa  : 

*  Complex  nonlinear  analysis  requiring  enormous 
amount  of  computer  time  due  to  the  iterative 
nature  of  the  problem. 

*  Solution*  to  very  complex  mechanical  devices 
consisting  of  many  links,  gears,  cams  and 
other  mechanical  elements. 


*  Handle  geo  me  r  r  i <:  .  deformation,  t»cr«-ss  a:.  *. 
frequency  constraints 

*  Requires  a  small  number  ot  function 
evaluations  1  iterations’  to  obtain  t  ri  •  primal 

design 

Supercomputer  needs: 

*  Combining  nonlinear  analysis  and  optimization 
involves  a  cons i derab 1 e  amount  of  computer 
time  for  each  solution. 

*  To  handle  objective  functions  which  are  highly 
nonlinear  functions  of  the  design  variables. 

Estimated  reduction  of  computer  time: 

100  to  I 


Estimated  reduction  of  computer  time: 

100  to  1 


2  .  E  f i  i  c  lent  Opt lmi rat  ion  Methods  for  Nonl inear 
Mech*n i sms  Des ign 

Objective:  Develop  efficient  optimization  methods 

with  the  following  characteristics  to  handle 
nonlinear  mechanism  design  problems: 

*  Minimize  an  objective  function  composed  of  a 
combination  of  rigid  body  and  deformation 
d i sp 1 acement • 

•  '§<“  geometric  and  cro  •  s-se  ct  ion*  l  size  design 
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Nonlinear  Analysis  and  Design 
of  Flexible  Mechanisms 

K.  D.  Willmert,  and 
M.  Sathyamoorthy 
Clarkson  University 

A  finite  element  approach  is  presented  here  for  the  nonlinear  vibrational  analysis  of 
planar  mechanisms.  The  analysis  takes  into  account  the  effects  of  material  and  geometric 
nonlinearities  on  the  dynamic  behavior.  The  geometric  nonlinearities  included  in  this 
-indy  are  due  to  stretching  of  the  neutral  axis  and  the  curvat ure-displacement  nonlinearity, 
both  caused  bv  large  deformations.  The  material  nonlinearity  is  due  to  a  nonlinear  stress- 
strain  relationship  of  hermile  polynomials  which  ensure  compatibility  of  curvative  between 
elements.  I  sing  a  variable  correlation  table,  a  global  system  of  nonlinear  equations  is 
derived  in  terms  of  the  global  unknowns  and  the  kinematics  of  the  mechanism.  A  harmonic 
series  t echnique  is  then  used  to  obtain  the  steady  state  solutions  to  this  system  of  nonlinear 
equations.  Numerical  results  are  presented  for  an  example  mechanism  and  the  effects  of 
the  nonlinearities  are  discussed. 

An  optimization  technique  which  is  applicable  to  problems  consisting  of  nonlinear 
objective  functions  and  constraints,  such  as  the  one  presented  here,  lias  also  been  dec  eloped 
are,  reported  I'lic  technique,  called  the  Causs  Nonlinear!)  Const  rained  Technique,  is 
i  llust  rated  w  it  h  example- . 
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THE  GAUSS  OPTIMIZATION  METHOD  FOR 
PROBLEMS  WITH  GENERAL  NONLINEAR  CONSTRAINTS 


T.  E.  Potter,  K.  D.  Willmert  and  M.  Sathyamoorthy 

Mechanical  and  Industrial  Engineering  Department 
Clarkson  University 
Potsdam,  NY  13676 


ABSTRACT 

A  new  algorithm  is  presented  Tor  -finding  the  minimum 
o-f  a  nonnegative  objective  -function  subject  to  general 
nonlinear  constraints.  This  algorithm,  based  o-f  Gauss' 
method  Tor  unconstrained  problems,  is  developed  as  as 
extension  to  the  Gauss  constrained  technique  -for  linear 
constraints.  The  derivation  of  the  algorithm,  using  a 
Lagrange  multiplier  approach,  is  based  on  the  Kuhn-Tucker 
conditions  so  that  when  the  iteration  process  terminates 
these  conditions  are  automatically  satisfied.  A  feasible 
design  is  maintained  throughout  the  iteration  process.  The 
solution  of  preliminary  examples  indicate  excellent  results 
in  terms  of  the  number  of  objective  function  evaluations 
required  by  the  algorithm  to  obtain  an  optimal  design. 


INTRODUCTION 

The  optimal  design  of  many  complex  structural  and  mechanical 
systems  is  hindered  by  the  large  computational  times  involved.  Most 
currently  available  optimization  techniques  require  a  large  number  of 
analyses  to  obtain  the  optimal  design.  For  small  problems,  or  ones 
in  which  the  analysis  is  simple,  these  methods  are  adequate;  however, 
for  large  problems,  or  where  a  time  consuming  anal ysi s  is  required, 
more  efficient  optimization  methods  are  needed.  The  goal  of  this 
research  was  to  develop  such  methods,  particularly  techniques 
applicable  to  mechanical  mechanism  design  where  the  members  are 
deforming  because  of  high  speed  motion  and  large  external  forces. 
Computat i onal  times  to  perform  a  single  analysis  are  enormous  for 
problems  of  this  type  involving  large  deformations  with  nonlinear 
material  char acter i st 1 cs.  Thus  the  goal  of  the  methods  developed  was 
to  reduce  the  number  of  analyses,  even  at  the  expense  of  increased 
: omputat l onal  effort  in  the  optimization  technique  itself,  i.e. 
additional  effort  in  finding  new  candidate  design  points  to  analyze. 
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'he  methods  developed  take  advantage  of  the  special 
;naracteristics  of  the  optimization  problem,  similar  to  the 
optimality  criterion  techniques.  This  greatly  improves  their 
itticiency.  For  most  mechanism  design  problems,  the  objective 
-unction  can  be  -formulated  as  a  sum  of  squared  quantities  such  as 
the  difference  between  the  desired  performance  and  the  actual 
performance  of  the  mechanism  at  specified  paints  during  its  motion. 
'hus  the  techniques  were  developed  specifically  to  handle  problems 
of  this  type,  although  the  methods  are  applicable  to  objective 
functions  which  are  general  sums  of  nonnegative  quantities,  such  as 
weight.  Many  mechanism  problems  have  constraints  which  are  only 
linear  functions  of  the  design  variables.  Thus  a  special  method  was 
developed  for  problems  of  this  type.  Other  problems  have 
constraints  which  are  linear  or  quadratic,  and  another  method  was 
developed  for  this  case.  Some  mechanism  design  problems  have  more 
general  nonlinear  constraints.  Methods  to  handle  these  cases  are 
currently  being  developed. 

All  of  the  techniques  developed  in  this  work  have  been  based 
on  Gauss’  method  CIO  which  is  applicable  to  problems  without 
constraints.  Wilde  C21  has  shown  this  method  to  be  particularly 
efficient  on  simple  mechanism  design  problems.  The  research 
presented  in  this  paper  has  extended  this  method  to  handle  various 
types  of  constraints  common  to  more  complex  mechanism  design 
probl ems. 


FORMULATION 

For  an  unconstrained  sum-of-squares  objective  function 


f  ( x  >  =  +  4* , 


( 1 ) 


wnere  <*>  is  a  vector  of  linear  or  nonlinear  functions  ♦  .  thru  ♦  in 

-*■  l  P 

:,  the  Gauss  method^for  calculating  the  next  iteration'of  the 
design  variables,  x  ,  given  a  current  design,  x  ,  is: 

K  1  K 


k+1 


-»  T 

J  <  x  )  J  ( x  ) 
k  k 


-1 


J  ( x  ,  )f(x,  )  , 
k  k 


(2) 


where 


J  (x)  =  V  ( x  >  .  .  .  V4>  (x) 


L  ' 


) 


=  v*T. 


(3) 


It  is  observed  that  only  first  derivatives  of  the  ♦  functions  are 
required  and  that  the  new  design  point  is  calculated  directly  from 
the  current  design  without  using  a  step  length  deter mi nat i on  with 
associated  one  dimensional  minimization.  Himmelblau  Cll  has  shown 
*:his  method  to  be  very  efficient  for  unconstrained  minimization 
jr  ob l ems . 


w: 


"his  technique  has  been  extended  to  handle  linear  inequality 
constraints  of  the  form 

—  -±T  ^ 

g  ( x  )  =  b  x  -  c  i  0 ,  1  =  1 . M  (  4 ) 

I  1  l 

as  well  as  equality  constraints 

h  ( x )  =  d x  -  e  =  0 ,  i  =  l . L .  ( 5  > 

i  l  i 

In  the  derivation  of  the^opt l mi z at i on  method,  the  4,  functions  are 
assumed  to  be  linear  in  x  of  the  form 

♦  =  JTx  u.  (&> 

where  J  is  a  constant  matrix.  However  the  resulting  technique  is 
applicable  to^problems  in  which  the  *.’s  are  general  nonlinear 
functions  of  x. 

At  iteration  k,  the  L  equality  constraints  and  any  of  the 
inequality  constraints  that  are  active  can  be  combined  and  written 
in  the  form 

BTx  -  C  =  0.  (7) 

If  at  the  next  iteration,  k  +  1,  the  variables  *k+1  are  at  the  optimum 
design,  then  the  Kuhn-Tucker  conditions  will  be  satisfied 

Vf  (x,  ,  )  BX  =  0  (8) 

K+  1 

BT?  4  -  C  =  0  (9) 

k  + 1 


X  >0  (10) 

where  X  is  the  vector  of  Lagrange  multipliers.  The  gradient  of  f  is 
given  by 

7f(x)  =  2J+(x>.  (ii: 

-k.  -X. 

Expanding  +(x)  in  a  Taylor  series  results  in 


♦  ( x ,  «  )  =  ♦ <x , )  + 

k  +  1  k 


JT[xk+1  -  *„]  (higher  order  terms)  (12) 


It  is  noted  £hat  the  higher  order  terms  are  equal  to  zero  if  *  is 
linear.  If  ♦  is  not  linear  then  these  terms  will  be  neglected  and 
the  expansion  is  only  approximate^  If  equation  (12)  is  substituted 
into  equation  (11),  evaluated  at  the  result  is 


=  2J  *(V  *  jT["k  +  i  ■  xk]  • 


This  may  be  substituted  into  the  first  ) uhn-Tucter  condition. 
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lave  the  -form 


1  -*■  -  -T- 

g  ( x  )  =  -xA  x  B  :<  -  C  £  0.  1  =  1  ,  .  .  -  ,  M. 

l  2  l  l  l 

It  at  iteration  1+1  there  are  r  active  constraints  (r  £  M> ,  the 
Kuhn-Tucker  conditions  will  be 


Vt(x  )  r  X  [a 
1  j  =  l 


k  +  1 


A  =  O 


J  J  J 


(IB) 


(19) 


and 


1  -*•  T  -  -  T  -*■ 

.Ax,  ,  +  B  x,  .  -  C  =0,  j  =  1 . r 

—k+ljk+l  j  k  +  1  j 


A  >  0 


(20) 


(21  ) 


where  the  summation  in  equation  (19)  and  the  j  subscript  in  equation 
(20)  refer  to  the  set  of  active  constraints  only. 


Using  a  derivation  similar  to  that  for  linear  constraints, 
substituting  the  expression  for  the  gradient  of  f,  equation  (13), 
^nto  the  first  Kuhn-Tucker  condition,  equation  (19),  and  solving  for 
produces 


k  ■+■  1 


“  ■ 
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T  V 

T  t  -*• 

2  J  J  -*■  )  A  X 

2JJ  x  -  Vf (x. )  -  )  B  X 

j  =  l  J  J 

k  k  “ .  J  j 

L  j-1  J 

(22) 


This  expression  for  x^  +  j  in  terms  of  A  is  now  substituted  into  the 
second  Kuhn-Tucker  condition,  equation  (20),  to  obtain: 
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T  V 
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2JJ  +  >  A  X 
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(23) 


These  r  nonlinear  equations  m^terms  of  the  new  unknowns,  X^  thru  X^_ , 
and  the  old  design  variables,  x  ,  are  solved  by  an  iterative  process 
for  the  values  of  *  thru  A  .  The  lambda  values  are  then  substituted 
into  equation  <22'  which  wiTl  /ield  new  values  for  the  design 


variables.  It  is  observed  that  the  matr  1  2JJ  is  the  matrix  of 
second  partial  derivatives,  G,  of  ^.he  objective  function  if  it  is 
quadratic.  Thus,  by  replacing  2JJ  in  equation  (22)  and  (23),  this 
technique  becomes  a  modification  of  the  second  order  method  rather 
than  the  Gauss  method. 

At  the  optimum  design  all  constraints  will  either  be  satisfied 
(less  than  zero)  or  active  (equal  to  zero)  and  each  active  constraint 
will  have  a  corresponding  lambda  whose  value  is  greater  than  or  equal 
to  zero.  If,  at  some  iteration  the  set  of  design  variables  yields  a 
violated  constraint,  then  obviously  the  optimum  point  has  not  been 
reached.  In  this  case,  the  newly  violated  constraint  will  be  added 
to  the  set  of  active  constraints  and  the  procedure  allowed  to 
continue.  If  at  some  iteration  the  set  of  design  variables  yields 
all  active  or  satisfied  constraints,  but  one  or  more  of  the  active 
constraints  has  a  corresponding  negative  lambda,  then  the  optimum 
design  has  also  not  been  reached.  The  negative  lambda  implies  that 
the  iteration  process  would  like  to  move  away  from  the  correspondi ng 
constraint  boundary  toward  the  feasible  region  where  the  constraint 
is  satisfied.  Thus,  the  constraint  is  dropped  from  the  set  of  active 
constraints  and  the  process  allowed  to  continue.  If  more  than  one 
negative  lambda  existed,  then  constraints  are  dropped  one  at  a  time 
starting  with  the  constraint  with  the  most  negative  lambda. 

A  constraint  is  added  to  the  set  of  active  constraints  if  it 
should  become  either  active  (equal^to  zero)  or  violated  (greater  than 
zero)  when  the  step  is  taken  from  x^  to  x^  .  In  the  case  where  a 
constraint  becomes  violated,  a  line  is  "drawn"  between  and  *k+1 
and  the  actual  step  is  taken  to  the  farthest  point  along  the  line  so 
that  no  constraints  are  violated.  In  effect,  this  procedure  is  the 
same  as  stepping  back  from  *k>1  toward  x  until  the  newly  violated 
constraint  is  just  active  (equal  to  zero).  The  constraint  is  then 
added  to  the  set  of  active  constraints  for  the  next  iteration. 

An  example  problem  with  quadratic  constraints  given  by  Boston, 
Willmert  and  Sathyamoor thy  C41  shows  this  method  to  be  very  efficient 
when  compared  to  the  generalized  reduced  gradient  method  (GRG) .  The 
problem  consisted  of  finding  the  optimal  design  of  a  four— bar 
mechanism  (minimizing  the  coupler  point  path  error  with  respect  to  a 
given  path)  subject  to  several  linear  constraints  on  link  length  and 
movability.  Additionally,  constraints  were  placed  on  the  crank  pin 
to  limit  its  location  to  the  intersection  of  two  circular  (quadratic) 
regions.  The  program  was  run  for  a  four  by  four  matrix  of  problems 
which  included  four  different  starting  points  and  four  different 
conditions  of  the  quadratic  constraints.  For  all  sixteen  runs,  the 
number  of  objective  function  evaluations  for  this  new  method  ranged 
from  8  to  32  (average  was  15),  while  the  GRG  method  required  from  303 
to  699  (average  was  502)  evaluations. 

The  interesting  information  here  is  that  the  solution  to  this 
quadr at l cal  1 y  constrained  four-bar  mechanism  used  no  more  objective 
^unction  evaluations  than  the  linearly  constrained  four-bar  mechanism 
e  ample  by  b'aradis  and  Willmert.  While  these  two  examples  are 
necessarily  different,  this  tendency  toward  requiring  similar  numbers 
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of  function  evaluations  for  different  classes  of  problems  is  very 
desirable.  The  net  result  is  that  we  now  have  an  optimization 
procedure  for  sum  of  squared  quantities  objective  functions  subject 
to  linear  and  quadratic  constraints  that  not  only  requires  relatively 
few  function  evaluations,  but  seems  to  be  constraint  order 
independent.  Now  the  need  is  to  determine  a  method  which  will  also 
work  for  higher  order  constraints. 


Boston,  et  al .  C4]  attempted  to  apply  this  method  to  higher 

order  problems,  but  met  with  mixed  results.  The  problems  encountered 
seemed  to  tie  in  with  the  higher  order  constraints  rather  than  with 
the  higher  order  objective  functions.  There  are  several  limitations 
implicit  in  the  algorithm  which  appear  to  be  the  source  of  the 
problems  encountered.  The  first  limitation  has  to  do  with  the 
application  of  the  constraints,  equation  (18),  to  the  first 
Kuhn-Tucker  condition,  equation  (19).  When  appr ox  1 mati ng  a  higher 
order  function  by  a  quadratic  Taylor  series  expansion  about  some 
point  x  ,  not  only  is  the  matrix  a  function  of  x^,  but  so  is  the 

B  vector  and  the  C.  scalar.  Thus  the  constraint  approx  1  mat i on , 
equation  (18),  shouid  be  written  as 


-  1-“T 

g  ( x  )  =  -x  A  ( x  .  )  x 

l  2  l  0 


[vv]  *  -  V; 


V 


£  0, 


=  1 


.  .  ,M 


(24) 


where 


■+•  ax  . 


As  x  approaches  (or  Ax  approaches  0)  then  this  approximation 
approaches  the  exact  value  of  the  constraint.  Thus  as  the  algorithm 
progresses  along  and  constraints  are  added  and  dropped,  the 
constraints  must  be  r eappr ox i mated  at  the  latest  design  to  keep  the^ 
step  size  small.  This  can  be  achieved  by  taking  the  new  values  of  x 
as  generated  by  equation  (22)  and  substituting  them  intg  the^ac^ual 
constraint  equations  to  get  improved  values  for  the  A^ (x^) ,  (x^) , 

and  (x  >  terms  in  equation  (24)  with  respect  to  the  current  design 
poi nt . 


The  second  limitation  involves  the  second  Kuhn-Tucker  condition, 
equation  (gO)  ,  which  is  used  to  obtain  the  equation  for  the  new 
values  of  X,  equation  (23).  This  is  simply  the  equation  for  the 
active  constraints.  In  the  original  formulation,  an  attempt  was  made 
at  obtaining  a  linear  approx i mat l on  in  X  for  this  constraint 
equation.  This  would  allow  equation  (23)  to  be  solved  explicitly  for 
X.  However,  failing  this  an  iterative  procedure  was  employed  to  find 
the  values  for  X.  Now  that  an  iterative  process  is  required  there  is 
no  advantage  in  keeping  a  quadratic  appr ox l mat i on  when  the  actual 
constraint  will  work  just  as  well.  Replacing  equation  (20)  with  the 
active  nonlinear  constraint  equations  will  remove  any  errors  due  to 
the  appr ox i mat l on  process. 

The  stepping  back  procedure  for  violated  constraints,  described 
above,  can  also  be  a  source  of  problems.  With  non-convex  programming 
problems  this  procedure  may  lead  to  a  situation  where  the  algorithm 
cannot  move  away  from  a  non -opt i mum  design.  Because  the  stepp l ng 


bact  procedure  assumed  a  straight  line  path  between  the  two  design 
points,  it  is  possible,  when  backing  out  of  a  newl /  violated 
constraint,  to  move  into  the  violated  region  of  the  constraint  that 
was  active  at  the  beginning  of  the  step.  The  procedure  would  then 
step  back  still  further  until  all  constraints  are  satisfied.  It  is 
possible  to  end  up  with  the  same  set  of  active  constraints  as  at  the 

start  of  the  iteration.  In  this  case  the  next  iteration  will  produce 

the  same  design,  which  may  be  non-optimal. 

Two  alternatives  are  readi 1 y  apparent  which  may  solve  this 
problem.  The  first  one  is  that  when  a  constraint  becomes  violated, 
repeat  the  step  but  include  the  newly  violated  constraint  in  the  set 

of  active  constraints.  The  second  alternative  is  to  move  to  the 

point  where  the  constraint  is  violated,  and  then  iterate  from  there 
without  stepping  back.  Of  course,  the  violated  constraint  is 
added  to  the  set  of  active  constraints.  Boston,  et.al.  C4D  looked 
into  this  second  alternative  to  some  extent.  They  reported  that  it 
did  not  always  work.  However,  it  is  not  clear  if  it  was  the  "no 
stepping  back"  that  was  the  cause  of  the  problems  or  if  the  second 
order  appr ox  1  mat i ons  to  the  constraints  contributed  to  the  problems. 


In  summary,  the  Gauss  nonlinearly  constrained  technique  is  very 
effective  at  solving  quadratical 1 y  constrained  problems.  No  major 
difficulties  appear  to  exist  which  would  preclude  it  from  solving 
problems  with  higher  order  constraints  once  the  modifications 
discussed  above  are  implemented.  This  method  with  the  proposed 
modi f ications  is  currently  the  leading  candidate  as  the  best  method 
for  solving  highly  nonlinear  mechanism  design  problems. 


RESULTS 


A  verification  of  the  effectiveness  of  the  Gauss  constrained 
method  applied  to  problems  with  quadratic  constraints  is  obtained  by 
solving  the  Rosen-Suzuki  test  problem  C5D: 


-k.  2  2  2  2 

minimize  Fix)  =  x +  x  2x^  +  x*T  -  5x  -  5x_  -  21x^  +  7x  . 


subject  to: 


g  t  (x  ) 


1 


+  x+  x,  +  x.  +  X,  -  X_ 


x  —  x  .  —  85:0 

3  4 


2 

y 

'  1 


2x~  +  x^  +  2x^  ~  x  ,  -  x  -1020 

2  3  4  1  4 


gT  <x  > 


=  r> 


2x 


+  *Z  +  2*1  "  *2  "  ;<  4  ~  5  S  0 


The  optimum  design  for  this  problem  is  at  x  =  [0,1,2,-lj 


Two  versions  o+  the  Gauss  nonlineariv  constrained  technique  and 
the  generalized  reduced  gradient  method,  identified  as  GRG,  were  used 
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tram  four  different  starting  points.  One  version  of  the  Gauss 
nonlinearly  constrained  technique,  identified  as  GNLC,  uses  the 
stepping  back  procedure  and  requires  a  feasible  starting  design  and 
will  always  maintain  a  feasible  design.  The  other  version, 
identified  an  GNLC. NS,  does  not  use  the  stepping  back  procedure  and 
has  no  requirement  on  the  feasibility  of  the  design  at  any  stage  of 
the  op 1 1  mi z at 1  on .  The  results  are  summarized  in  Table  I.  It  can 
easily  be  seen  that  the  Gauss  nonlinearly  constrained  technique  is 
much  more  efficient  with  respect  to  number  of  function  evaluations 
than  the  generalized  reduced  gradient  method. 


ALGORITHM 

STARTING 
DESIGN,  k 

NUMBER  OF 
ITERATIONS 

FUNCTION 

EVALUATIONS 

GNLC. NS 

[0,0,  0,0] 

2 

~7 

GC.NS 

[1,  1, 1,  1] 

T 

GC .  NS 

[2, 2, 2, 2] 

2 

3 

GNLC. NS 

[0,0,7§,0] 

*T 

4 

GNLC 

[0,0,75,0] 

3 

4 

GNLC 

[0,0, 0,0] 

5 

6 

GNLC 

[1,1.1,  1] 

5 

6 

GRG 

[0,0,75,0] 

9 

83 

GRG 

[0,0, 0,0] 

1  1 

106 

GRG 

[1,1,1,  1] 

1 1 

133 

GRG 

[2, 2, 2,2] 

1 1 

144 

Table  I:  Comparison  of  Algorithms 


CONCLUSIONS 

The  optimization  techniques  developed  in  this  research  as 
extensions  of  the  Gauss  method  to  handle  various  types  of 
constraints  are  effective  approaches  toVeducing  the  number  of 
analyses  required  to  obtain  an  optimal  design.  As  a  result,  the 
computational  time  for  1 arge  problems  should  be  reduced 
s i gn if tc an  1 1 y . 
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ABSTRACT 

A  finite  eleaent  approach  la  preaantad  In  chla  paper  for  the  nonlinear  vi¬ 
brational  analysis  of  planar  aechanlsaa.  The  analyala  takes  Into  account  the 
effecta  of  aaterlal  and  geosMtrlc  nonllnearltlea  on  the  dynaalc  behavior.  The 
geometric  nonllnearltlea  Included  In  thla  study  are  due  to  stretching  of  the 
neutral  axis  and  the  curvature-dlsplaceaent  nonlinearity,  both  caused  by  large 
dcfornatlone .  The  aaterlal  nonlinearity  la  due  to  a  nonlinear  atreaa-straln  re¬ 
lationship  of  the  Haaberg-Oagood  type.  The  analysis  presented  here  aakea  use  of 
heraite  polynoalals  which  ensure  coapatibillty  of  curvature  between  eleaents. 
Using  a  variable  correlation  table,  a  global  syetea  of  nonlinear  equations  are 
derived  In  terms  of  the  global  unknowns  and  the  klneautlcs  of  the  aerhanlsa .  A 
harmonic  series  technique  Is  then  used  to  obtain  the  steady  stats  solutions  to 
thla  system  of  nonlinear  equations.  Numerical  results  are  presented  for  an  ex¬ 
ample  mechanism  and  Che  effecta  of  the  nonllnearltlea  are  discussed. 

INTRODUCTION 

The  Importance  of  flexibility  of  linkages  on  the  performance  of  high  speed 
alnlatas-auiss  aechanleaa  la  well  recognized.  A  considerable  aaount  of  research 
has  been  done  In  this  area  In  the  last  two  decades.  While  It  Is  desirable  to 
develop  enalytlcal  and  numerical  procedures  that  enable  the  design  of  rigid  link 
aechanlsaa  and  robots  to  perform  a  given  function  with  specified  reliability.  It 
la  also  Important  to  evaluata  Che  effects  of  flexibility  of  elastic  aeabers  on 
on  their  performance.  It  la  known  that  a  mechanism  designed  for  operation  at 
low  speeds  may  not  perform  satisfactorily  at  high  speeds  due  to  Che  effects  of 
large  Inertia  forces  and  resulting  elastic  deformations.  Thus  It  becoaes 
necessary  to  Include  In  the  dynaalc  analysis  of  aechanlsaa,  not  only  the  effect 
of  the  rigid  body  action,  but  also  the  flexibility  of  the  linkages. 

Moat  of  the  previous  Investigations  In  the  area  of  elastic  analysis  of 
aechanleaa  have  been  carried  out  within  the  framework  of  Che  linear  theory  (1- 
171.  However,  Vlacoal  and  Ayre  (181  used  a  Galerkln-type  nonlinear  analysis 
procedure  to  study  the  vibrations  of  a  slider-crank  aechanisa.  A  later  work  by 
Sadler  and  Sandor  [191  used  the  lumped  paraaeter  approach  to  a  nonlinear  dynaalc 
aodel  of  an  elastic  llnksge.  The  aechanisa  analyzed  In  this  paper  was  a  general 
four-bar  linkage  and  the  analytical  aodel  included  the  response  coupling  associ¬ 
ated  with  both  the  transmission  of  forces  at  the  pin  joints  and  the  dependence 
of  the  undeformed  notion  of  a  link  on  the  elastic  motion  of  other  links.  A 
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finite  element  jnaLvsis,  with  the  aLd  of  the  piecewise  linear  method  of  Martin, 
was  used  by  Sevak  and  McLarnan  (201  to  carry  out  the  nonlinear  analysis  of  a 
mechanism.  Further  nonlinear  work  dealing  with  the  vibrations  of  elastic  mecha¬ 
nisms  are  reported  in  References  (21-23).  In  a  recent  investigation,  Thompson 
and  Sung  (24)  used  a  variational  formulation  for  the  nonlinear  finite  element 
analysis  of  planar  mechanisms  considering  geometric  non  1 inea r i t ies .  Some  experi¬ 
mental  results  were  also  presented. 

This  paper  is  concerned  with  the  nonlinear  vibrational  analysis  of  general 
planar  mechanisms.  A  finite  element  method  is  used  which  includes  the  effects  of 
both  geometric  and  material  nonlinearities.  The  geometric  non l inea ri t ies  in¬ 
cluded  in  this  study  are  due  to  stretching  of  the  neutral  axis  with  partially 
constrained  ends  and  a  general  curvature-displacement  relationship,  both  caused 
bv  large  deformations.  The  material  nonlinearity  is  of  the  Ramberg-Osgood  type 
with  three  parameters  to  represent  the  nonlinear  s c re ss - st ra in  relationship  (25- 
27].  Additional  effects  considered  are  transverse  shear  and  rotatory  inertia 
and  changes  in  cross-section  due  to  realistically  proportioned  members.  The 
governing  nonlinear  differential  equations  are  derived  for  each  element  in  terms 
of  the  axial  and  transverse  deformations,  rotations,  curvatures,  and  shear  defor¬ 
mation  angles.  These  equations  are  then  assembled  with  the  aid  of  a  variable 
correlation  table  and  the  resulting  global  system  of  equations  is  solved  using  an 
iterative  technique  based  on  a  harmonic  series  solution  procedure. 

FINITE  ELEMENT  FORMULATION 


A  finite  element  method  is  presented  below  for  the  nonlinear  analysis  of  a 
general  closed  looped  mechanism.  The  mechanism  can  be  composed  of  various  com¬ 
binations  of  simple  four  bar  chains,  frame  elements,  sliders  moving  on  fixed  ref¬ 
erences,  or  sliders  moving  on  rotating  links.  Each  link  is  divided  into  one  or 
more  elements  with  each  element  having  the  local  coordinate  system  as  shown  in 
Figure  1.  If  a  slider  is  present,  the  masses  and  M2  are  located  at  ends  1  and 

2,  as  shown.  The  length  of  the  element  A  is  constant  except  for  links  with 
9lider9  moving  along  them. 

The  displacement  vector  of  any  point  (a)  on  the  element's  neutral  axis  is 
given  by: 

S  *  (X^  cosy  +  siny  +  x  +•  u)i  +  (Y^  co9y  -  X^  siny  +  v)j  (1) 

where  X[  and  Yj  are  the  coordinates  of  end  1  of  the  element  given  by  the  rigid 
body  motion.  The  coordinate  x  is  measured  along  the  element's  neutral  axis  from 
1  to  2  and  y  is  the  angle  between  the  rigid  body  position  and  the  X-axis.  The 
axial  and  transverse  displacements  of  point  (a)  from  the  rigid  body  position  are 
given  bv  u  and  w,  respectively.  This  equation  takes  into  account  both  the  rigid 
body  motion  and  the  elastic  displacements  and  defines  the  position  of  any  point 
along  the  neutral  axis. 

Differentiating  Equation  1  with  respect  to  time  yields  the  velocity  of  any 
point  (a).  The  unit  vectors  i  and  j  move  with  coordinate  system  and  vary  with 
time.  The  angular  velocity  of  any  differential  line  segment  on  the  neutral  axis 
of  the  element  is  given  by: 


where  y,c  i9  the  derivative  of  y  with  respect  to  time  and  w,xt  is  the  derivative 
of  the  transverse  displacement  with  respect  to  the  local  coordinate  x  and  time  t. 
The  kinetic  energy  due  to  rotation  of  the  element  is  given  by: 


h  r  1 

f  f^|oAx| S,tl ^  +  plz(y,t  +  w,xt)2idx  dy 


-  7  ^1 1 s • t ! 2  +  I  O' 

2  x  *0  2  *-  x*A 

where  0  13  the  mass  density  and  Ax  and  Iz  are  the  cross  sectional  area  and  moment 
of  inertia  of  the  element  respectively.  The  term  oIz(v,t  +  w,xC>v2  represents 
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the  effect  of  rotatory  inert ic.  The  kinetic  energy  due  to  beam  bending  associ¬ 
ated  with  transverse  shear  [28!  is: 


■  I  2 1 ,  £  }dx  dy 


where  i  is  a  measure  of  the  transverse  shear  angle. 
The  strain  energy  for  the  element  is  given  by: 


"  *  .  vol  dE  dvo1  *  ,’vol  7  Txv  Yxv  dvo1  (5) 

o  *- 

where  ",  ,  'xv,  and  vxv  are  the  normal  stress,  normal  s  rain,  shear  stress  and 

shear  strain,  respectively.  For  a  nonlinear  material  of  the  Ramberg-Osgood  type 
[25-4.71,  the  relationship  between  stress  and  strain  is: 


where  A  corresponds  to  Young's  modulus  E,  and  Bem  represents  the  nonlinear  term. 
A,  3  and  ra  are  constants  for  the  particular  material  being  considered.  The  above 
relationship,  Equation  6,  is  valid  only  for  positive  strains.  If  the  strain  is 
negative,  the  following  expression  is  used: 

7  3  A*:  +  B  (-c>m  if  r  0  (7) 


The  change  in  sign  of  the  nonlinear  term  results  in  the  same  overall  effect  on 
the  s t re ss - s t ra in  relationship  as  for  positive  strain,  i.e.  either  hardening  or 
softening  depending  on  the  values  of  B  and  ra. 

Using  the  Ramberg-Osgood  re lat lonship ,  the  following  expression  for  the 
strain  energy  is  obtained  for  positive  strain  z  ^  0: 


h 

•  2  -.1 


‘1  A  f2  '  i^I  ®  c^1'dx  dY  +  'h  ’A  \  Gxy  Y xy  dx  dy 


where  x  is  the  axial  coordinate,  v  is  the  transverse  coordinate,  and  Gxy  is  the 
shear  modulus.  When  r  0  the  equation  is: 


- — r  B  (-c)»*Mdx  dv  + 

2  nr*-l  -  h 


'h  2  Gxy  Y xy  dx  dy 


The  nonlinear  expression  for  the  curvature  R  of  a  planar  static  beam  under¬ 
going  large  deformations  is: 


(1  +  w  2)3/2 


Thu9  the  strain  is 


w,  2)3/2 


Combining  the  geometric  non  1  inea r i t ie s  due  to  stretching  of  the  neutral  axis  and 
the  curvature-displacement  nonl inearitv,  results  in  the  expressions  for  normal 
and  shear  strain: 
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(1  +  Db  w,x>3/2 


Substituting  these  expressions  for  strain  into  the  strain  energy  Equations 
(8)  and  (9)  produces,  for  e  >_  0  or  e  <0: 


,  ?  ,  A  ] 

!  . h  '  b  { -  A  ( u  .  x  + 

—  o 
2 

D>  1 

-  —  B  (±  u.x  ±  i 

m+t  2 


*,2)3/2 


-  Gxy  (w,x  +  a)2}dx  dy 


(14,15) 


For  equation  (15)  negative  sign  Is  applicable  in  all  the  term*  containing  +  signs 
and  la  valid  for  e  <  0.  In  Equations  (12),  (14)  and  (15),  D^,  and  D„  are 
tracing  constants  representing  Che  effects  of  material  non  1 Ineari t y ,  geometric 
nonlinearity  due  to  curvature,  and  geometric  nonlinearity  due  to  stretching  of 
neutral  axis,  respectively.  Each  tracing  constant  is  equated  to  unity  when  that 
particular  nonlinearity  is  being  considered  and  Is  equated  to  zero  when  it  is 
not . 

In  order  to  represent  realistically  proportioned  members,  changes  In  cross 
section  are  Included.  Each  element  la  divided  into  sections  of  varying  lengths 
with  constant  area.  The  integrations  Involved  In  the  element  equations  are 
carried  out  In  a  piecewise  fashion  with  the  area  in  each  section  taken  as  a  con¬ 
stant.  This  procedure  provides  a  reasonable  approximation  of  variable  cross 
sectional  umbers  without  having  to  resort  to  large  numbers  of  elements. 

The  Lagranglan  function  L  is  defined  as 
S  «k 

L  .  I  l  (K.E.g  +  K.E.g  -  5.E.  )lk  (16) 

k-1  1=1 

where  Nk  is  the  total  number  of  elements  In  the  kth  link  and  S  Is  the  number  of 
links  in  the  mechanism.  Substituting  Equations  (3),  (4)  and  114)  Into  Equation 
(16),  the  Lagrangian  .  can  be  expressed  in  terms  of  the  displacements  u  and  w, 

Che  shear  angle  a,  and  the  rigid  body  toot  ion. 

Hermite  polynomials  are  used  to  approximate  u,  w,  and  a  in  order  to  satisfy 
the  boundary  conditions  of  various  types  of  mechanisms  easily  and  to  ensure  in¬ 
terelement  compatibility.  The  axial  deformation  u  is  approximated  by  a  linear 
shape  function  given  by 

u(x,t)  -  U^t)  Ni(x)  +  U2(t)  N2(x)  (17) 

Similarly,  fifth  degree  polynomial  shape  functions  are  used  to  approximate  the 
transverse  deformation  v: 

v(x,t)  ■  Wj(t)  Hi}(x)  +  9 1  ( t )  H2i(x)  +  *i(t) 

+  W2(t)  Hi2(x)  +  02(t)  H22(x)  +  m2(t)  H32(x)  (18) 

The  shear  angle  a  is  also  approximated  by  a  fifth  degree  polynomial  in  order  to 
make  it  compatible  with  the  transverse  displacement  w.  Therefore  a  is  assumed 

to  be : 


a(x,t)  ■  a^(t)  Hjk|(x)  +  *i(t)  H2i(x)  +  X^(t)  H3j(x) 

+  a2(t)  Hi2(x)  +  a>2(t)  H22(x)  +  X2(t)  H32(x) 
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where  *.  and  X  are  Che  first  and  second  derivatives  of  a,  respectively. 
The  Hemite  polynomials  are  given  by: 


N  l  (  x ) 

•  l  -  e  ;  N2(x) 

■  e 

Hn(x) 

-  1  -  I0e3  +  15e4  - 

6e3  ; 

Ht \ ( x ) 

■  A(e  -  be3  +  Se^  - 

3e5)  ; 

H  3 1  ( x  ) 

•  -  3e3  +  3e^* 

-  e5>/2 

(20) 

H i 2 ( * )  "  10e^  -  I5e^  +  6e^ 

"  ^(-4e3  +  "  3e^) 

;  H32(x)  *  42(e3  -  2 tu  +  e5)/2  (21) 


where,  e  *  x/A. 

A  transformation  of  coordinates  is  now  introduced  to  change  from  the  moving 
coordinate  system  associated  with  the  elements  to  global  coordinates.  Only  U]_, 
Wj,  l?2  and  W2  need  to  be  transformed.  The  other  coordinates  are  angles  or  de¬ 
rivatives  of  angles  which  are  not  directional  on  the  X,  Y  coordinate  system  used. 
The  transformations  are: 


l!l  •  m  cos$i  -  W|  sin$i  ;  U2  *  U2  cos*2  “  w2  8in$2 

W j  *  U 1  sini^  +  Wi  ;  W2*U28in$2+^2C08^2  (22) 

For  pin  connections  the  transformation  angles  and  $2  are  set  equal  to  -y  (the 
rigid  bodv  angle)  which  transforms  the  coordinates  back  to  the  global  coordi¬ 
nates.  For  sliders  moving  on  rotating  links  the  transformation  becomes  more  in¬ 
volved.  In  this  case,  the  deformation  of  the  driver  link  must  be  transformed  to 
correspond  to  the  axial  and  transverse  deformation  of  the  output  link  (14). 
Substituting  the  expressions  from  Equation  (22)  into  Equations  (17)  and 
13),  the  global  coordinates  for  the  system  are  then: 

q  “  !l'j  W|  ^  nj  aj  K/i  U2  W2  §2  ai2  32  ^2  ^2^  (23) 


The  Lagrangian  function  is  then  written  in  terms  of  the  transformed  element  co¬ 
ordinates.  Differentiating  the  Lagrangian  with  respect  to  the  element  coordi¬ 
nates,  the  following  element  equations  are  obtained: 


_d_ 

dt 


(4^-) 

*q .  t 


3L  _ 


0 


(24) 


In  differentiating  the  expressions  for  the  kinetic  and  strain  energies  in  Equa¬ 
tions  (3),  (4)  and  (14),  it  must  be  remembered  that  A,  which  is  the  upper  limit 
of  integration,  is  a  function  of  time.  The  operations  carried  out  In  Equation 
<24)  results  in  a  system  of  nonlinear  element  differential  equations.  Assem¬ 
bling  the  element  matrices  for  the  particular  mechanism  being  solved  results  in 
the  global  system  of  equations: 

M  Q,tt  ♦  C  Q,c  +  (Ke  +  Kn)  Q  -  F(t)  (25) 

The  M,  C,  Ke  and  matrices  are  all  functions  of  time.  The  C  matrix  results 
from  the  kinetic  energy  of  the  system.  No  damping  was  Included  in  the  formula¬ 
tion  of  the  problem.  The  C  Q,c  term  was  found  to  be  small  and  thus  was  ignored 
in  the  analysis.  The  matrix  Ke  is  the  linear  portion  of  the  total  stiffness 
matrix.  It  is  a  function  of  rigid  body  motion,  but  is  not  a  function  of  the 
deformations  Q.  The  matrix  ,  however,  is  the  nonlinear  portion  of  the  stiff¬ 
ness  matrix.  It  is  a  function  of  the  deformations  Equation  (25)  is  thus  a 

nonlinear  system  of  differential  equations. 

The  derivation  of  the  finite  element  Equation  (25)  Is  based  on  the  assump¬ 
tion  of  positive  strains  e.  If  the  strain  is  negative  a  similar  derivation  is 
possible,  based  on  Equation  (15)  for  the  strain  energy  rather  than  Equation 
14).  The  only  difference  in  the  resulting  Equation  (25)  is  in  the  stiffness 
matrix.  Wherever  an  ®  ^  occurs  in  Che  original  formulation,  it  becomes  (-e)®~^ 
for  negative  strains.  All  other  negative  signs  resulting  from  the  introduction 
of  -*  cancels  out  in  the  d i f f erent iat ions  required.  Thu9  in  order  to  handle  both 
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positive  and  negative  strains  the  terms  involving  in  the  stiffness  matrix 

were  replaced  by  -e  ra'^. 

In  order  to  solve  the  nonlinear  svstern  of  Equation  (25),  an  iterative  ap¬ 
proach  was  used.  First  the  equations  were  solved  using  the  linear  terms  onlv, 
i.e.  Che  matrix  was  ignored.  This  was  accomplished  bv  setting  all  of  the 
tracing  constants  Db  ,  Dg  and  equal  to  zero.  The  solution  Q  for  the  linear 
equations  was  then  used  to  determine  values  for  the  nonlinear  stiffness  raatrLx 
.  Equation  (25)  was  Chen  solved  again  for  new  values  of  Q,  and  the  process  re¬ 
peated.  Experience  showed  chat  this  procedure  converged  in  from  3  to  5  itera¬ 
tions.  To  solve  Equation  (25)  for  Q( t ) ,  for  particular  Kn ,  a  harmonic  series 
solution  method  was  used  similar  to  that  of  Bahgat  and  Willmert  !l4).  This  ap¬ 
proach  overcomes  problems  with  stability,  due  to  the  time  varying  nature  of  the 
matrices,  that  sometimes  result  from  an  eigenvalue  technique.  The  steadv  state 
solution  is  obtained  without  adding  artificial  damping.  The  solution,  without 
the  C  Q,c  term,  is  given  by: 

N 

Q(t)  *  (K  •  n“-  (A.  cos  n-t  +  sin  n-t)  (26) 

n*o  n 

where  w  is  the  input  crank  speed,  and  and  Bn  are  solutions  to  the  linear  equa- 


Fftu)  *  L  co9  nuttv  +  m  Bn  sin  nut'*  ;  k  *  0 ,  l ,  .  .  .  ,  2.N- 1  (27) 

*  nmo  ^  *  n*o  n  K 

where  B0  is  set  equal  to  zero.  The  values  of  t-K  are  the  times  at  2N  equal  time 
increments  per  revolution  of  the  input  crank  given  by: 

for  k  =  0,1 . 2N-1  (28) 

Computational  experience  indicates  chat  a  fairly  accurate  solution  is  obtained 
using  only  a  few  terms  in  Equation  (26).  As  the  number  of  terms  increases  the 
components  of  the  matrix  (K  -  n^  <u^M)  grow  and  thus  Che  inverse  (K  -  n^  ^M)"* 
becomes  small.  The  summation  can  therefore  be  truncated  to  reduce  computational 
t  ime . 

The  stress  in  the  links  is  calculated  by  evaluating  the  strains  from  Equa¬ 
tion  (12).  The  stress  can  then  be  determined  at  any  point  in  an  element  using 
Equations  (6)  and  (7).  To  find  the  maximum  stress  in  an  element  the  maximum 
strain  must  be  found.  Setting  the  first  derivative  of  Equation  (12)  to  zero  and 
solving  the  resulting  expression,  the  position  of  the  maximum  strain  is  deter¬ 
mined.  Once  the  location  is  known,  the  maximum  strain  and  stress  can  be  evalu¬ 
ated  . 

The  above  formulation  is  based  on  the  use  of  the  shear  angle  a,  which  is  ap¬ 
propriate  particularly  for  short  members.  For  long  slender  links  this  quantity 
is  not  required.  The  elimination  of  a  reduces  the  size  of  the  problem  consider¬ 
ably  sLnce  the  nodal  deformations  aj,  X^,  12,  *2  <*nd  ^  would  no  longer  be 

present.  For  long  slender  members: 


Using  this  expression,  the  equation  for  strain  energy  (14)  for  positive  shear  e 
reduces  to: 


,2  ,A  rl  l 

Jhi  b  '2 A  <Ui* +  2 


(l  +  Db  *.2)3/2 


(U- 


(1  +  Db  w, 2)3/2 


)dx  dy 


A  similar  expression  exists  for  negative  strain.  The  kinetic  energy  also  changes 
if  a  is  not  present.  The  energy  associated  with  transverse  shear,  Eouat ion  (4), 
is  eliminated  and  thus  Equation  (3)  represents  the  total  kinetic  energy  of  the 
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element.  Using  a  procedure  similar  to  Che  method  outlined  above,  vibrational 
equations  of  che  same  form  as  Equation  (25)  can  be  obtained,  but  they  will  be 
smaller  in  size.  However,  nonlinear  terms  still  exist  in  the  stiffness  matrix 
due  to  che  material  and  geometric  non l inea r 1 c les .  The  method  of  solution  is  thus 
identical  to  that  outlined  earlier. 

EXAMPLE  PROBLEMS 

The  following  example  is  presented  to  illustrate  the  method  of  solution. 

The  nonlinearities  due  to  neutral  axis  stretching,  and  complex  curvature-dis¬ 
placement  and  stress* strain  relationships  are  all  considered.  A  four-bar  link¬ 
age,  as  shown  in  Figure  2,  is  used  as  the  example  with  all  of  the  members  flexi¬ 

ble  and  made  of  the  same  material.  The  data  for  the  mechanism  is:  Length  of 
input  crank  (AB)  •  5.0  in.  Length  of  coupler  (BC)  ■  11. 0  in,  Length  of  rocker 
(CD)  ■  10.5  in,  Fixed  distance  (AD)  •  10. 0  in,  Cross  section  of  links  a  rec¬ 
tangular,  Height  of  rectangle  *1.0  in,  and  Width  of  rectangle  *  0.25  in.  The 
position  of  the  input  crank  is  zero  degrees  at  t  ■  0  and  the  direction  of  rota¬ 
tion  is  counterclockwise.  The  mechanism  la  divided  Into  three  elements  with  each 
link  in  the  mechanism  taken  as  an  element.  The  boundary  conditions  are  that  only 
moment  and  shear  terms  exist  for  the  input  crank's  driven  end  (A).  For  the  pin 
connections  between  links,  there  are  deformations,  rotations  and  shear  terms,  and 
for  the  rocker's  fixed  point  there  are  only  rotation  and  shear  terms. 

First,  the  deformations  in  the  mechanism  were  determined  with  the  shear 
angle  a  present .  In  this  case  the  crank  link  was  rotated  at  100  rad/sec.  The 

material  properties,  approximating  aluminum,  were  as  follows:  A  -  10.87xl06 

lb/in^,  B  *  0.8387x10^'  lb/in^,  m  *  3.0,  and  Mast  density  *  0.0002536 
lb-sec^/in^.  Three  separate  procedures  were  used  to  obtain  numerical  results. 
First  che  problem  was  solved  using  the  linear  analysis  method  of  Bahgat  and 
Willmert  (14] ,  with  E  ■  A.  Next  the  method  of  this  paper  was  used  with  the  trac¬ 
ing  constants  equal  to  zero.  Thus  a  linear  analysis  was  obtained.  Finally  the 
method  was  applied  with  all  tracing  constants  equal  to  one,  i.e.  a  full  nonlinear 
analysis.  A  representative  deformation  Uj  as  a  function  of  crank  position  ia 
shown  in  Figure  3.  This  Is  the  horizontal  deformation  of  the  free  end  of  the 
crank  link.  As  can  be  seen,  che  three  curves  are  very  similar.  The  effect  of 
the  shear  angle  a  is  to  increase  the  deformation  slightly.  For  thia  slow  speed 
the  linear  and  nonlinear  analyses  were  almoat  identical. 

The  same  problem  was  solved  again  at  a  higher  speed  of  200  rad/sec.  The  re¬ 
sulting  deformation  Uj  is  9hovn  in  Figure  4.  As  can  be  seen,  high  frequency 
oscillations  started  to  appear,  with  greater  separation  between  the  three  analy¬ 
ses.  At  even  higher  speed  these  oscillations  became  more  predominant  to  che 
point  of  instabilities  in  the  motion  at  very  high  speeds. 

The  revised  form  of  the  analysis  equations  was  considered  next,  i.e.  the 
form  without  the  snear  angle  a.  Here  a  crank  speed  of  150  rad/sec  was  used.  A 
comparison  was  made  of  the  effects  of  the  various  non l inear 1 t les  on  the  deforma¬ 
tions  and  stresses  as  compared  to  the  linear  analysis.  Figures  5  and  6  show  a 
comparison  of  the  linear  and  nonlinear  deformations  U2  (the  horizontal  deforma¬ 
tion  of  the  free  end  of  the  output  link)  caused  separately  by  geometric  and 
material  nonlinearities.  Figures  7  and  8  show  the  maximum  stresses  in  the  con¬ 
necting  link  of  the  mechanism.  As  expected,  the  material  nonlinearity  of  the 
Ramberg-Osgood  type  results  in  deformations  which  are  greater  in  magnitude  than 
those  obtained  using  a  linear  elastic  model.  The  maximum  stress  decreased  due  to 
the  presence  of  che  term  Bem  substracted  from  the  linear  stress  expression. 

The  geometric  nonlinearities  considered,  namely  curvature  displacement  and 
stretching  of  the  neutral  axis,  both  due  to  large  deformations  produced  mixed 
results  with  deformations  reduced  at  9ome  points  and  increased  at  other  points. 
The  effect  of  the  geometric  nonlinearities  would  be  expected  to  produce  a  stiff¬ 
ening  of  the  members  [28]  of  the  mechanism  and  thus  produce  smaller  deformations. 
The  increased  deformations  in  this  case  might  be  due  to  the  fact  that  the  deform¬ 
ations  are  in  relationship  to  the  entire  mechanism  and  not  just  to  an  individual 
beam  element. 

CONCLUSIONS 

The  nonlinear  analysis  procedure,  using  a  finite  element  technique,  is  an 
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effective  method  of  calculating  the  steady  state  deformations  and  stresses  In  a 
mechanism.  Significant  differences  can  occur  between  the  linear  and  nonlinear 
approaches.  This  was  particularly  true  for  the  stresses  in  the  example  consid¬ 
ered  in  this  work.  Research  Is  still  needed  on  the  overall  effect  of  the  shear 
angle  a,  and  a  more  complete  picture  of  the  nonlinear  terms  In  Che  analysis  would 
be  of  value.  Additional  nonlinear  effect  should  also  be  investigated,  such  as 
the  effect  on  the  translations  of  one  link  due  to  large  deformations  of  the  other 
1  Inks . 
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Presented  in  this  paper  is  a  new  optimization  technique,  called  the  Gauss  nonlinearlv  constrained 
technique,  which  is  applicable  to  desutn  problems  with  nonlinear  objective  functions  and  constraints  The 
technique  is  an  extension  of  a  previously  developed  method  for  linear  constraints,  referred  to  as  the 
Gauss  constrained  technique  Both  of  these  techniques,  based  on  the  Gauss  unconstrained  method,  have 
been  developed  so  that  the  Kuhn-Tucker  conditions  are  automatically  satisfied  when  the  procedure 
terminates. 


1.  Introduction 

The  optimal  design  of  many  structures  and  mechanical  mechanisms  involves  one  or  more 
complex  and  time-consuming  analyses  at  each  iteration  of  the  optimization.  This  is  particularly 
critical  if  a  large-deformation  nonlinear  analysis  is  required.  In  these  cases  especially,  it  is 
important  that  the  optimization  method  requires  very  few  analyses,  even  at  the  expense  of 
significantly  increasing  the  amount  of  calculations  by  the  optimization  technique  itself.  For 
example,  in  a  recent  work.  DeRubes  and  Willmert  [1]  applied  the  relatively  efficient  generalized 
reduced  gradient  (GRG)  technique  of  Lasdon  et  al.  [2]  to  mechanism  design  for  path 
generation  and  rigid-body  guidance.  The  mechanism  links  were  considered  flexible,  and  thus  a 
quasi-static  (linear)  finite  element  analysis  was  used  to  obtain  deformations  and  stresses.  The 
GRG  required  as  many  as  14.000  mechanism  analyses  to  obtain  the  optimal  design.  Comput¬ 
ation  time  approached  twenty  hours  on  an  IBM  4341  mainframe  computer.  If  a  nonlinear 
analysis  had  been  used,  the  corresponding  times  would  have  been  considerably  higher. 

To  reduce  the  number  of  analyses.  Paradis  and  Willmert  [3]  developed  a  new  direct  method 
for  efficient  design  of  mechanisms.  Gauss'  method,  which  Wilde  [4]  concluded  to  be  very 
efficient  for  unconstrained  mechanism  design,  was  modified  to  handle  linear  constraints.  The 
resulting  technique,  referred  to  as  the  Gauss  constrained  technique,  was  highly  efficient  and 
required  very  few  objective  function  evaluations  to  obtain  an  optimal  design.  Their  method  has 
been  extended,  in  this  paper,  to  handle  nonlinear  constraints. 
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2.  Development  of  the  method 

The  optimization  problem  consists  of  minimizing: 

F( x)  =  if' if  .  (  1 1 

where  ip  is  a  vector  of  general  functions  of  the  variables  x  Many  optimal  design  problems  have 
objective  functions  of  this  form.  For  the  derivation  of  the  optimization  technique,  the  ip  are 
approximated  by  linear  functions  of  x  (which  makes  F  quadratic)  of  the  form: 

ip  —  y 1  x  -  c ,  i  - ) 

where  J  and  c  are  a  constant  matrix  and  vector,  respectively.  However,  the  technique,  once 
derived,  will  be  aplied  to  more  general  cases,  where  the  <p  are  highly  nonlinear  functions  of  x. 
The  constraints  are  approximated  by  a  general  quadratic  of  the  form: 

g  (x)  =  ix'A  x  ~  b'x  -  </  d  .  /  =  1 . k  .  ( 3  ) 


The  gradient  of  the  objective  function  (  1  )  is  given  by 

'FF(x)  =  IJip  .  (4) 

which  is  exact,  whether  ip  is  linear  or  not.  as  long  as  J  is  the  matrix  of  lirst  partial  derivatives  of 
ip.  The  matrix  of  second  partial  derivatives  of  the  objective  function  is 

a  =  :jj'.  (5) 

which  is  exact  only  when  ip  is  linear  At  any  given  iteration  we  assume  there  are  /  active 
constraints  (ordered  such  that  /  =  1 .  2.  .  .  / ).  where  /  -s  A\  Thus,  the  Kuhn-  Fucker  conditions 


V/dx)  -  (  A  x  -  b  (A.  - 


;  A  x  *-  b  )  A  =d 


( x )  =  ;  x  A  x  b  x  -  ii  -il.  j  =  1 . 


A^(i.  /  = 


Similar  to  the  development  of  Paradis  and  Willmert  |3).  the  new  method  is  generated  such 
that  a  single  iteration  yields  the  optimal  design  for  a  quadratic  objective  function,  assuming  that 
the  constraints  active  at  the  optimal  design  are  also  active  at  the  starting  point.  If  F{x)  is 
quadratic,  then  the  following  equation  is  valid  for  anv  two  points  x  .  and  x  : 

V/(x  .  )  FF(x  l  -  (/ix  x  I  (d| 

Solving  (  d  |  for  C Fix  .  I.  assuming  r  is  the  optimal  design,  and  substituting  into  ( results  in 
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the  iterative  expression  (after  solving  for  x^  . .  ): 

x..|-|G*i  .-A  A  )  Gx,  -  ^  b  A  -  VF{xt  )  (ID) 

;  L  .  .  j 

If  the  objective  function  has  the  special  form  of  ( 1 ).  then  the  matrix  C,  in  the  iterative  equation 
would  be  replaced  by  (5).  i.e.  G  =  2 JJ\  Otherwise,  the  method  could  still  be  applied  as  long  as 
the  matrix  of  second  partial  derivatives  G.  or  an  approximation  to  it.  were  known 

In  order  to  use  the  iterative  equation  (  10).  the  unknown  vector  of  Lagrange  multipliers  A 
must  be  determined.  This  is  accomplished  by  substituting  the  expression  for  x  . ( 10).  into 
each  of  the  active  constraints  of  (  7 )  resulting  in  a  system  of  l  nonlinear  equations  in  /  unknowns 
A  of  the  form: 


A)  =  jj(G'I  A  A.) 

i  -  J 

Gx  -  j]  b  A,  -VF(x.  ) 

L  -  I 

x  j(G  -  S  A  A.) 

•  -  I 

i  /- 

Gxt  —  ^  b,\  -  VF(xt ) 

*  b  \  (  G  J-  A ,  A  ) 

- 1  r 

Gx„  -  £  M,  -^F(x. 

L  -  1 

This  system  of  equations  in  A  could  be  solved  using  several  different  methods.  The  approach 
used  in  this  research  was  as  follows.  At  each  iteration,  initially  the  Newton-Raphson  technique 
was  applied.  If  the  method  did  not  converge,  which  may  be  caused  by  the  fact  that  the  equations 
had  no  solution  or  that  the  method  simply  was  not  able  to  locate  one.  then  another  approach 
was  used.  In  this  case,  an  objective  function  H  was  formed  which  was  the  sum  of  the  squares  of 
the  active  constraints  g, .  This  unconstrained  function  was  minimized  with  respect  to  A.  In  this 
work  Powell's  method  was  used,  but  any  available  technique  could  be  applied.  The  optimal 
values  of  A  (whether  H  is  zero  or  not  for  these  values  of  A)  were  then  used  in  the  iterative 
equation  (  10). 

As  in  the  original  Gauss  constrained  technique  (as  well  as  the  gradient  projection  method), 
the  criterion  for  dropping  a  constraint  from  the  set  of  active  constraints  is  the  sign  of  A,.  At  any 
iteration  the  constraint  corresponding  to  the  most  negative  A  is  dropped  from  the  set  of  active 
constraints  and  a  new  vector  of  Langrange  multipliers  is  determined.  This  procedure  is  repeated 
until  all  A,  2*  0.  If  this  results  in  no  active  constraints,  then  the  optimization  technique  continues 
by  setting  all  A,  =0  in  (10).  In  this  case,  the  iterative  expression  can  be  shown  to  reduce  to 
Gauss'  method  for  unconstrained  optimization: 

x  . ,  =  xt  —  ( yy 1  >  lj<p  (12) 

At  each  iteration  a  decision  must  be  made  as  to  whether  or  not  a  new  constraint  should  be  added 
to  the  set  of  active  constraints  If  the  direction  of  minimization  v,  .  is  defined  as  the  direction 
from  x  tor.  .  then  a  constraint  is  added  onlv  if  a  step  in  the  s,  . ,  direction  will  not  satisfy  the 
constraint 
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[wo  different  methods  were  used  to  handle  \iolated  constraints  The  tirst  approach  was  to  do 
nothing  special  it  a  constraint  became  violated  as  a  result  ot  an  iteration  based  on  (KM  Thus, 
whether  a  constraint  is  only  active,  re.  g  =  M  or  violated,  i.e  g  •  <•.  it  is  treated  the  same  in  the 
optimization  technique  Ihe  basis  behind  this  approach  was  the  assumption  that  eventually  the 
method  will  vitistv  all  constraints,  since  these  are  equations  (111  which  the  method  attempts  to 
satistv  at  each  iteration  This  approach  allows  the  optimization  tit  start  at  any  design  whether  it 
satisfies  the  constraints  or  not. 

The  other  method  ot  handling  violated  constraints  was  to  start  at  a  point  which  satisfied  all 
constraints  Then,  it  a  step  is  taken  using  (  KM  such  that  g,(x  .  )  Ml.  a  new  design  x  _  was 
found  along  the  line  connecting  x  and  x  .  such  that  g.txj  ..)  =  ().  This  was  accomplished  hv 
stepping  back  from  a  violated  constraint  using  the  approximate  equation 

(x  .  -  x  KMx  i 

x  -  x -  (13) 

g,lx,  .  I  CM*  l 

Ihe  design  x  .  was  then  used  in  the  next  iteration. 

Although  the  technique  was  developed  assuming  linear  functions  v?(x).  i  e.  F (x)  quadratic, 
the  method  is  applicable  to  problems  where  the  y?  are  general  functions  of  x.  An  especially 
important  characteristic  for  mechanism  design  is  that  the  technique  requires  a  total  of  only  one 
more  objective  function  evaluation  than  iterations  to  obtain  an  optimal  design,  unlike  the  GRG 
technique  of  Lasdon  et  al.  [2j.  where  a  step-length  determination  is  required  at  each  iteration 


3.  Example  problems 

Several  examples  were  considered  in  this  work,  two  ot  which  are  presented  here. 

The  first  problem,  which  is  a  modified  Rosen-Suzuki  test  problem,  consists  of  a  quadratic- 
objective  function  in  eight  variables  with  seven  quadratic  inequality  constraints  as  follows: 

Min  F(x)  =  xT  -  •  x]  -  :  i  -  3xj  -  2x:  -  ;  x."  ~  x  •  -  :  x ; 

-  2.x.  -  A i  -  5  v .  -  v<  +  "  x„  -  2  x  x,  . 

S.t.  g . (X)  =  i .x 7  \x\  -  2x\  -  \x~-  -  x;  -  2x.  4 1,  -  5 v -  -  KMM)  , 

g . ( x )  =  x;  -  r.  -  v:  -  x\  ~  xl  *•*  *;  ~  t .  -  G  -  r .  x. 

*  v  i  .  -  i  i ,  -  i  -  x\  fs  -  ii. 

g  ( x  i  =  i  -  2 x  -  i  -  2  x  .  -  x  ~  2  x .  *  x  •  2  i  .  i  x .  x  v ,  13*  1 1  . 

g  lx)  -  x  -  2  i  -  ;  i  -  4 » ;  -  5  x .  -  h x  -  i  x  x  \  .  s  •  1 1 

g(  X I  -  X  -  I  .  X  1  .  I  I  •  • 1  . 
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Cjx)  -  2 v '  2v,  -  2aA  2 v 4  ~  3x;  -  3x;  -  3v:  -  3a; 

-  2a  -2a,  2a. .  -  2x  -  2a-  -  2a,  -  1(1===  0  , 

c-ixi  -  4 1 -  4v;  -  2a;  -  2a;  -  3  a.  -  5a,  -  5  a,  -*•  5v-  -  20^0  . 

fables  1  and  2  show  the  results  of  the  optimization  from  two  different  starting  points  Three 
different  techniques  were  applied.  The  first  two  were  variations  of  the  Gauss  nonlinearly 
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constrained  method  The  hrst  one.  labeled  GNLC'.NS  in  the  tables,  ignored  violated  con¬ 
straints.  while  the  second,  labeled  GNLC.  treated  violated  constraints  using  1  13).  Both  of  these 
methods  were  compared  with  the  generalized  reduced  gradient  technique,  labeled  GRG 

Table  1  shows  the  results  from  a  starting  point,  where  no  constraints  were  active  or  violated 
All  three  methods  produced  the  same  optimal  design.  The  number  of  iterations  NI  and  number 
of  function  evaluations  NT  for  the  GNLC  methods  were  considerably  less  than  for  the  GRG 
Table  2  shows  the  results  from  a  starting  point  in  which  three  constraints  were  violated,  while 
one  was  active  For  this  starting  point  only  the  GNLC. NS  technique  was  applicable.  Again,  the 
number  of  iterations  and  function  evaluations  were  less  than  for  the  GRG. 

The  second  example  consist .,  of  the  design  of  a  four-bar  mechanism  such  that  the  coupler 
point  i  XX.  YY)  of  the  mechanism  closely  generated  the  curve  defined  by  the  eight  points  as 
shown  in  Fig.  1  The  objective  function  to  be  minimized  is  the  sum  of  the  distances  (squared) 
between  the  desired  curve  and  the  actual  curve  generated  bv  the  mechanism: 


w  here 


Fix)  =  ^  (  XX  XG  r  *  (  YY  -  YG  )'  . 

XX  =  v,  cos  y  -  v<  cos  1 1)  -  v,„  .  YY  =  ,v.  sin  y  *  ,v<  cos  <i/  *  x ,4 

<!>  -  -  v-  -  t?  -  y  . 

i .  sint  y  -  v„ ) 


(-■  tan 


T)  =  cos 


L  x.  -  x  2  cos(  y  -  )  J 

G  *  X]  -  vt  -  x\  -  2.r,.t;  cos(  y  -.rjt 


-V.  v, 


i :  -  Zx  x .  cos(  y,  -  ,t„ )] 1 
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v  -  -  Ay 

and  t\Cj.  VG  I  are  the  coordinates  of  the  desired  curve. 

I  he  design  variables  x. . v,,.  are  link  lengths,  orientation  of  some  of  the  links  as  well  as 

the  location  of  the  crank  pin.  Note,  the  objective  function  is  a  highlv  nonlinear  function  of  the 
design  variables.  This  is  the  same  problem  that  Paradis  and  Willmert  [3]  solved,  however,  two 
nonlinear  constraints,  g,,  and  g.4.  have  been  introduced  here  to  constrain  the  location  of  the 
crank  pm  to  two  user-defined  circular  regions.  Thus,  our  design  problem  contains  ten  variables, 
r  .  .  v  ,.  and  fourteen  inequality  constraints  g  (jc)«l).  /  =  1 . 14.  The  constraints  are: 


tr  ( x )  =  1 

-  i .  ^  n  . 

g:(x)  =  1  -  .v,  -s  0  . 

g..(X)  —  v; 

-  V,  s:  II  . 

g4(x)  =  x:  -  x,  -c  0  . 

c.(  X  1  =  V. 

-  v,  -  U  . 

J?«( x)  =  x,  *  x ,  -  x.  -  x 4  «  0  . 

e  - 1  x )  =  i . 

*  l  -  V  .  v  1 1  . 

y\(x)  =  x,  +  v4  -  x,  -  x,  0  . 

C.IX)  =  V 

3i ) 1 1  , 

g|n(x)  =  x,  -  3U  s;  0  , 

g  IX)—  V 

.  ■  3i  i  <  i )  . 

II 

1 

/A 

o 

e  :  1  X  )  =  i . 

v -  CX, )"  -  (  r.,  -  CY, )' 

g,4(x)  =  (xU)  -  CX;):  +  (xy  -  CY:): 

-/?;•£  o.  -/?;*;<). 
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Table  3 


Results  of  mechanism  design  problem 


Starting 

point 

GNLC 

GRG 

Status 

NT 

1 1 

feasible 

NF 

12 

389 

15.0 

.V. 

30  0 

30  0 

5.0 

.V. 

8.345 

8  345 

15.0 

•V. 

25  821 

25.821 

5.0 

V. 

12.817 

12.  sr 

15.0 

.V, 

16.026 

16.026 

-0.007 

-C 

-1063 

-1  063 

0.412 

AY 

1  259 

I 

1.379 

AY 

1.118 

1  118 

6.0 

v. 

4,0 

4  0 

5.0 

A'.,, 

'.0 

■’  ll 

220.29 

F 

1  '8349 

1 .  ’8348 

— 

<7*.  1 

9.13.14 

9.13.14 

'Angle  variables  (radians) 


where  CX,,  CY,.  /?,.  CX:,  CY,.  and  R ,  were  4.  4.  3.  7.  7.  and  3,  respectively. 

The  resulting  optimal  designs  using  the  GNLC  and  GRG  methods  are  shown  in  Table  3.  The 
optimal  mechanism  is  shown  in  Fig.  2.  Again  the  Gauss  nonlinearly  constrained  method 
required  significantly  fewer  objective  function  evaluations. 

4.  Conclusions 

The  Gauss  nonlinearly  constrained  method  is  an  effective  technique  of  solving  nonlinear 
design  problems.  It  is  particularly  efficient  for  cases  in  which  the  evaluation  of  the  objective 
function  is  very  time-consuming.  Although  in  some  instances  considerable  calculations  must  be 
done  per  iteration,  the  amount  required  is  still  insignificant  compared  to  that  required  to  do  just 
one  analysis  of  a  complex  mechanical  or  structural  system. 

The  method  has  been  applied  to  a  variety  of  problems  consisting  of  objective  functions  of 
various  complexities.  In  all  cases  it  has  worked  well.  However,  it  has  received  only  limited 
application  to  problems  involving  constraints  that  are  more  than  quadratically  nonlinear.  It 
appears  that  further  research  is  required  for  such  cases. 
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Rumberg-Osgood  tvpe  is  employed  for  the  nonlinear  stress-strain  behavior. 
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element  analysis,  with  high  order  heraite  polynomials  which  ensure  compati¬ 
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ferential  equations  are  solved  by  means  of  a  harmonic  series  technique  to 
obtain  steady  9tate  solutions.  The  effects  of  the  non  1 i nea r i t i e s  are  dis¬ 
cussed  by  means  of  an  example  problem. 


Aft 


i» 


No n  —  'M inrar  111 i li r  a  1 t  o us 
S>labUth|,  anb  Dijumiiirs 
S’lnix'litrrs  aub  .'Hr <  ban  is 


23-25  March  1987 

A hs Ira c I s 


Sponsored  hij 
Army  Research  OfJ'i('r 

1 1  os  Led  ()}j 

Virginia  I ^olylechn w  Insliluh 
and  .S’  la!  e  // ni  i  ua'si  I  // 


Effect  of  Geometric  and  Material 
Nonlinearities  on  Vibration  of  Planar 

Mechanisms* 


Iv  B.  kear  111  M.  Sat liyamoorl  It y 

K.  1).  Will mort 

Depart  iiiriil  of  Meeliauu  al  and  Imlu. dri.il  l  .ii^iueering 
Clarkson  I'niversify 
I’otsdam,  New  York  KU)7f> 


Abstrac  t 

A  nonlinear  vibrational  analysis  for  planar  mechanisms  with  ma¬ 
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Figure  1:  Global  Displacements  For  Element  i 
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Figure  2:  Deflected  Element. 


•  A'i  is  the  A'  coordinate  of  end  1  of  the  reference. 


•  \\  is  the  V  coordinate  of  end  1  of  the  reference. 


o  is  the  angle  between  the  reference  and  X  axis. 


•  n  is  the  axial  displacement  of  a  point  from  its  reference  position. 

•  u:  is  the  transverse  displacement  of  a  point  from  its  reference  position 


•  r is  a  unit  vector  in  x  direction. 


/"is  a  unit  vector  in  y  direction. 
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LSODE  (Livermore  Solver  for  Ordinary  Differential  Equations)  can  be  use 
■olvc  a  system  of  equations  of  the  form 
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Consider  a  straight  beam  element  with  its  elastic  axis  parallel  to  the  x-axis  of  the  fixed 
global  system  of  coordinate  G.  Under  the  action  of  loads  w  hich  are  distributed  along  the  beam 
element,  it  deforms.  This  deformation  is  described  by  the  displacement  of  each  point  of  the  elastic 


axis.  The  components  of  these  displacements  u,  v  and  w  are  in  the  direction  of  &xG,  CyG  and  e/G,  re¬ 
spectively  (see  Figure  1 ).  In  addition,  each  cross  section  is  also  rotated  by  angle  0  about  the 
elastic  axis. 

Due  to  this  deformation,  the  triad  cxj_[,  evj_j  and  ez^  at  the  point  B  of  the  undeformed  axis  is 
transformed  into  a  new  rotated  orthogonal  triad  of  unite  vectors  Cy^,  ezP)  and  the  position 
vector  of  the  point  B,  an  clastic  axis,  after  deformation  becomes 

RBD  =  ,XG  +  Ql-^xG  +  ( YG  +  Qn^yG  +  (ZG  +  Q3)6zG 
+  (xp  +  u)exH  +  v'byH  +  w5z.H 

while  the  position  vector  RpD  of  the  point  P  off  the  elastic  axis  after  deformation  reduces  to 

Rp  =  (XG  +  Q 1  )<fxG  +  ( YG  +  Q2^yG  +  (ZG  +  Q3)(5zG 
+  (Xp  +  u)exH  +  veyH  +  we7H 

+  yp  gyp  +  /p&zD  +  :UxD 

w  here  a.  is  the  warping  displacement  of  the  cross  section  which  is  considered  an  extension  of  the 
classical  St.  Vcnant  torsion  to  a  more  general  case. 

To  this  point  the  elastic  deformatioas  of  the  element  arc  defined  for  the  general  case  where 
the  deformatioas  are  due  to  the  action  of  forces  which  arc  distributed  along  the  beam  element.  In 
the  case  of  mechanisms,  the  element  is  imposed  to  an  additional  angles  of  deformation  yx,  Yy  and  yz, 
due  to  the  elastic  deformations  of  the  other  elements,  in  the  directions  x^j,  y^  and  z.pj,  respective¬ 
ly. 

Finally  ,  the  position  vector  of  the  point  P  after  deformation  with  respect  to  the  undeformed 
■o  s tem  of  coordinates  w  ill  reduce  to 

Rp  =  L<Xr  +  Qai>  +  ,*H  +  uT|TrH 
+  '  xp  TDr  TrllU-  °H 

where 

xr  =  uxG  yg  zg  j 
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Qal  ~  tQl  $2  Q3J 

Xj-j  =  LXh  0  Oj 


-T 

U  =  lu  v  wj 

*p  —  J  p  ^p  J 

Trpj  is  the  elastic  transformation  matrix  and  is  given  by 

TrH  =  TrD  TDH  where 
TpH  >s  a  transformation  matrix  of  the  body. 

The  velocity  of  the  point  P  after  deformation  at  any  instant  of  time  is  given  by 

RpD  =  LrxT  +  Qal)TGH  +  (^T  +  §al)TGH 

+  ( x^  +  +  0  TrH  +  xp  f  rHj  eH 

+L(XT  +  Qa])TGH  +  (XH  +  €T)TrH  +  xpTDrTrPi)ji;H 
Using  finite  element  expression  to  rewrite  the  previous  equation  and  isolating  the  velocity  terms, 
the  expression  can  be  written  in  the  following  form 


Rp  Lax  atp  aqe  J 


,qe 


Also,  if  the  axial  deformations  in  the  direction  exj_j.  evj_j  and*z^  arc  transformed  to  the  global  fixed 
coordinate,  the  final  expression  for  the  velocity  reduces  to 

*^p  ~  Lax  a<p  ^  J 

ai  Kinetic  Energy 

If  we  consider  the  previous  beam  clement  with  an  additional  two  concentrated  masses  <mcon] , 
mc()n->)  at  nodes  1  and  nod2  respectively,  then  the  kinetic  energy  of  the  element  becomes: 

KH  =  1  nf>  RpRpdvol  .  1/2  R,  R, 

’•l  *  1/2  "W 

where  p  is  the  mass  density  of  the  clement,  mcon ]  mcon->  arc  the  concentrated  masses  of  the  nodes  1 

•  * 

arid  2,  respectively  and  Rj  R->  are  the  velocities  of  the  nodes  1  and  2. 


3 


The  kinetic  energy  ot  the  element  can  be  rearranged  in  the  following  form 


K=  l/2^T0qjM 


b)  Strain  Energy 

The  strain  energy  of  the  element  including  material  nonlinearity  is 


SE  =  J  [(1/2  Eexx) 

VM. 


2.  _i_  ,B  m+1; 

m+1  xx 


1-  1/2  G  +  e*z>ldvol 

Lagrangian  Expression  is  given  by 


L  =  KE  -  SE. 


c)  Examples: 

1 1  Rotating  cantilever  beam  (see  Figures  2-7) 


Length  =  20.0  in 


radius  =  0.2  in 


E  =  10.H7E6  lb/in“ 


p  =  2.536  r-4  lb. s/in 
operating  speeds 


to  =  I  (X).  rad/s,  3(X)  rad/s.  5(X)  rad/s 


Deformation  using  the  formula  given  by  'Stephen.’ 


500  rad/s 


6.3799  E-4 


5.7419E-3 


1.5949E-2 


1.3914E-4 


1.2523E-3 


3.4787E-3 


1 .3905  E-4 


-1.2513E-3 


-3.47648E-3 


2 1  Four-Bar  Linkage  (see  Figures  8-2 1 1 
Height  of  the  beam  =  1 .0  in 


Width  of  the  beam  =0.25  in 


E  =  10.87  E6  lb/in“ 


B  =  8  387  E10  lb/in2 


"•r  V  ~T  *V*  ~J>  V*  "LF 


Coordinates,  u;  =100.  rad/sec. 
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Figure  6  Axial  Displacements  in  Global  Coordinates,  u j  =500.  rad/ 
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re  13  Maximum  Absolute  Stress  of  the  element  2.  ^  =  100.  rad/sec 


18 


2/2 


-A198  644 


UNCLASSIFIED 


NONLINEAR  ANALVSIS  AND  OPTIMAL  DESIGN  OF  DVNAKIC 
MECHANICAL  SVSTEHS  FOR  SPACECRAFT  APPLICATIONS) 
CLARKSON  UNIV  POTSDAM  NV  K  D  MILLMERT  ET  AL  SEP  87 
AF0SR-TR-87-2888  AF0SR-84-8876  F/G  22/2 


NL 


ill  |.0  *-  it 


m  112.0 


L25  I  U  mi  1.6 


MICROCOPY  RESOLUTION  TEST  CHAR! 


.V*>1 


"V *> 


